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''Zr^The  thermochemical  sublimation  response  of  ATJ-S  graphite  in  both  low 
and  high  mass  transfer  convective  environments  was  studied  both  experi¬ 
mentally  and  analytically.  ATO-S  graphite  models  were  ablation  tested  in 
the  NASA  Ames  Advanced  Entry  Heating  Simulator  (AEHS)  test  facility  at  >0.1 
and  91.0  atmosphere  impact  pressure  test  conditions.  Ablation  tests  were 
run  under  both  subsonic  and  supersonic  flow  conditions  at  both  pressures. 
High  quality  recession  rate  and  surface  temperature  from  the  ablation  tests 
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were  correlated  with  various  thermochemical  ablation  predictions.  Results 
from  the  low  subsonic  flow  conditions  (low  mass  transfer  rates)  exhibit  good 
agreement  with  the  JANAF  thermochemical  model  for  carbon  vapor  species  (Ci  - 
C5).  Results  from  these  tests  under  supersonic  flow  conditions  (high  mass 
transfer)  show  the  data  to  exhibit  slight  kinetic  sublimation  effects  com¬ 
pared  to  JANAF  equilibrium  sublimation  predictions.  Extrapolation  of  the 
inferred  kinetic  sublimation  effects  to  the  high  pressure  regime  of  interest 
(10  -  300  atmospheres)  shows  the  uncertainty  in  both  surface  temperature  and 
recession  rates  to  be  less  than  1  percent.  The  conclusion  is  made  that  the 
JANAF  equilibrium  sublimation  model  is  most  appropriate  for  graphite  thermo¬ 
chemical  ablation  predictions  under  high  pressure  reentry  conditions 
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This  Final  report  documents  the  results  of  research  performed  by  the  Aerotherm  Division  of 
Acurex  Corporation,  Momtain  View,  California,  for  the  Air  Force  Materials  Laboratory,  Systems  Sup¬ 
port  Division,  Space  and  Missiles  Branch,  under  Contract  F3361 5-74-C-5D94.  The  AFML  Project  Engineer 
was  initially  Capt.  G.  V.  Jumper,  Jr.  (AFML/MXS)  and  this  responsibility  was  subsequently  assumed  by 
Mr.  G.  G.  Ormbrek  (AFML/MXS).  The  Aerotherm  Project  Manager  was  Mr.  C.  A.  Powars.  This  research 
was  conducted  during  the  period  April,  1974  through  January,  1976.  This  report  was  submitted  to 
AFML  in  January  1976. 

The  experimental  portion  of  this  program  was  accomplished  using  NASA  Ames  test  facilities, 
equipment,  and  personnel,  through  a  cooperative  arrangement  between  the  Air  Force  Materials  Lab¬ 
oratory  and  NASA  Ames  Research  Center.  The  NASA  Ames  Advanced  Entry  Heating  Simulator  facility  was 
used  to  conduct  the  tests,  and  certain  required  modifications  to  this  facility  were  provided  by  NASA 
Ames  and  Aerotherm  under  AFML  sponsorship.  The  cooperation  of  Mr.  Howard  Larson  (Branch  Chief)  and 
Dr.  Philip  Nachtsheim  (Assistant  Branch  Chief)  of  the' Thermal  Protection  Branch  in  making  available 
the  NASA  facilities,  equipment,  and  personnel  to  support  this  program  is  gratefully  acknowledged. 
Valuable  technical  assistance  was  provided  by  Mr.  John  Lundell.  Mr.  David  Peterson  and  Mr.  Frank 
Nichols  provided  assistance  relative  to  test  equipment  utilization  and  modification.  Finally,  the 
services  of  Mr.  Jack  Hagan  and  Mr.  Harry  Simmon  in  operating  the  test  facility  and  associated  equip¬ 
ment  is  most  gratefully  acknowledged. 
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PRSOEDIIG  PACAj^ANKwNOT  'FILMED 


StCTlON  I 
INTRODUCTION 

The  purpose  of  this  study  was  to  examine  the  sublimation*  response  of  graphitic  materials  in 
the  high  convective  mass  transfer  rate  regime,  i.e.,  the  regime  of  primary  interest  relative  to  the 
application  of  these  materials  to  reentry  vehicle  thermal  protection  systems.  The  specific  objec¬ 
tives  were  to  (a)  design  and  demonstrate  an  experimental  procedure  for  isolating  the  effects  of  sub¬ 
limation  kinetics  in  the  high  convective  mass  transfer  rate  regime,  (b)  analyze  the  resultant  data 
to  determine  the  importance  of  sublimation  kinetics  in  controlling  the  rate  of  graphite  ablation  in 
this  regime,  and  (c)  analyze  these  experimental  data  in  conjunction  with  the  results  of  other  recent 
investigations  to  better  establish  the  correct  equilibrium  carbon  vapor  pressure  law. 

Important  uncertainties  currently  exist  regarding  the  sublimation  of  graphitic  materials  due 
primarily  to  uncertainties  in  the  carbon  vapor  equilibrium  thermochemical  properties  (i.e.,  the  car¬ 
bon  equation  of  state)  and  the  quantitative  importance  of  sublimation  kinetics  relative  to  boundary 
layer  diffusion  in  controlling  the  material  ablation  rate.  Uncertainties  also  exist  regarding  carbon 
melting  and  particulate  mass  loss  from  graphitic  materials  during  ablation.  These  uncertainties  are 
reviewed  in  Section  2.  In  this  section,  the  relevant  findings  of  prior  investigations  of  carbon  sub¬ 
limation  are  reviewed  and  compared.  Also,  current  mathematical  modeling  of  graphitic  material  abla¬ 
tion  is  reviewed  (both  with  and  without  sublimation  kinetics).  This  serves  to  establish  the  require¬ 
ments  for  experimentally  investigating  sublimation  kinetic  effects  in  the  high  convective  mass  transfer 
regime. 

The  experimental  approach  used  to  study  graphite  sublimation  kinetic  effects  is  described  in 
Section  3.  Since  this  experimental  approach  was  both  unique  and  complex,  details  regarding  the  test 
facility,  instrumentation,  and  test  procedures  are  carefully  documented  in  this  section.  Section  3 
also  includes  a  detailed  summary  of  the  resultant  test  data.  An  analysis  of  ;:he  experimental  data  is 
presented  in  Section  4.  Consideration  is  given  to  the  implications  of  the  data  relative  to  both  the 
carbon  vapor  pressure  law  and  the  effects  of  sublimation  kinetics. 

•Throughout  this  report,  the  term  "sublimation"  is  used  (rather  than  "vaporization"  which  is  often 
used  in  the  carbon  literature)  to  emphasize  that  it  is  the  solid  to  gaseous  phase  change  which  is 
being  studied. 
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SECTION  2 


CARBON  SUBLIMATION  UNCERTAINTIES 

The  hyperthermal  ablative  response  of  carbon  and  graphitic  materials  continues  to  be  contro¬ 
versial  due  principally  to  uncertainties  in  (1)  the  carbon  equation  of  state,  (2)  the  chemical  and/ 
or  diffusional  rate  processes  controlling  sublimation,  and  (3)  the  significance  of  particulate  mass 
loss  from  high  temperature  graphite  surfaces.  This  section  will  present  an  overview  of  work  in  the 
field  of  high  temperature  carbon  thermochemical  response  from  an  engineering  perspective  Topics 
reviewed  will  include  (1)  carbon  sublimation  thermochemistry  (Section  2.1),  (2)  carbon  melting  and 
solid  polymorphic  phases  (Section  2.2),  (3)  carbon  sublimation  rate  uncertainties  (Section  2.3),  and 
(4)  particulate  mass  loss  from  polycrystalline  graphites  and  carbon/carbon  materials  (Section  2.4). 

2.1  CARBON  SUBLIMATION  THERMOCHEMISTRY 

Extensive  work  in  the  area  of  equilibrium  carbon  sublimation  has  been  conducted  within  the  past 
two  and  a  half  decades  (Reference  1-20).  The  purpose  of  these  studies  has  been  to  identify  the  signif¬ 
icant  carbon  vapor  species  and  to  evaluate  their  thermochemical  properties  in  the  hyperthermal  regime 
of  interest  (i.e.,  T  >5400°R).  The  procedure  generally  used  to  study  equilibrium  carbon  sublimation 
is  mass  spectrometric  analysis  of  the  vapor  species  effusing  from  a  Knudsen  cell.  Some  experimental 
programs  however,  have  studied  sublimation  as  free  vaporization  directly  from  the  surface  of  the  specimen. 

The  most  recently  completed  program  utilizing  the  Knudsen  cell  approach  was  that  of  Milne,  et  al . , 
(Reference  14)  under  AFML  sponsorship.  Milne's  apparatus  was  unique  in  that  the  Knudsen  cell  was  oper¬ 
ated  at  temperatures  approaching  6000°R.  Their  measurements  of  ion  species  ratios  at  1 7  eV  showed  good 
agreement  with  the  results  of  Drowart,  et  al . ,  (Reference  1)  for  species  (^-Cj  and  although  their 
values  for  C4  were  somewhat  lower.  Milne's  results  do  therefore  essentially  substantiate  the  carbon 
vapor  thermochemistry  reported  by  JANAF  (Reference  21).  These  results  in  turn  tend  to  refute  the  the¬ 
oretical  carbon  equations  of  state  of  Dolton,  et  al.,  (Reference  22)  and  Kratsch,  et  al.,  (Reference  23) 
recommended  several  years  earlier. 

Wachi  and  Gilmartin  (Reference  15)  also  studied  carbon  sublimation  using  the  technique  of  free 
vaporization.  These  researchers  studied  the  free  vaporization  of  ATJ,  ZTA,  TU-6,  and  pyrolytic  graphites. 
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Relative  ion  intensities,  apparent  rates  of  vaporization,  and  activation  energies  of  vaporization 
were  measured  for  species  -f.^  in  the  temperature  range  5000“ R  -  5400’R .  Some  differences  in  the 
iree  vaporization  ol  polycrystalline  and  pyrolytic  graphites  were  noted.  These  differences  were 
apparent  at  temperatures  a!  jve  5200^  where  the  polycrystalline  graphites  would  exhibit  rapid  in¬ 
creases  in  their  vaporization  rates  with  an  attendant  ejection  of  crystallites.  (Note  that  Milne's 
results  did  not  exhibit  this  behavior.)  Wachi  and  Gilmartin  attributed  this  behavior  to  the  two 
phase  nature  of  polycrystalline  graphites  although  they  could  not  provide  a  precise  explanation. 

Interestingly,  this  observation  of  increased  mass  loss  with  a  substantial  particulate  component 
has  been  observed  by  others  (References  19,  24,  25)  and  relates  to  the  uncertainty  of  carbon  melting 
which  is  reviewed  in  Section  2.2. 

Wachi  and  Gilmartin  conclude  from  their  results  that  for  conventional  graphites,  the  sublination 
of  carbon  species  from  the  crystalline  phase  is  the  predominant  process,  and  at  temperatures  above 
53Q0°R  •  go°R  the  vaporization  of  species  from  the  amorphous  carbon  binder  is  the  predominant  pro¬ 
cess.  In  addition,  Wachi  and  Gilmartin  observed  a  time  dependency  of  the  relative  ion  intensities 
of  carbon  species  at  a  given  surface  temperature,  which  they  attribute  to  changes  in  surface  area  and 
surface  morphology.  Their  findings  relate  to  the  subject  of  carbon  melting  in  Section  2.2  and  specif¬ 
ically  the  results  of  Whittaker,  et  al . ,  (Reference  19). 

In  summary,  the  results  of  their  research  substantiate  the  JANAF  thermochemical  models  for 
vapor  species  Cj-Cj..  Thus,  the  recent  carbon  sublimation  studies  whether  using  a  Knudsen  cell  approach 
or  free  vaporization  support  the  JANAF  thermochemical  model  for  the  vapor  species  C^Cg  (Reference  21). 

The  predominant  species  Is  with  C1  and  both  about  a  factor  of  5  below  in  concentration.  Spe¬ 
cies  and  Cg  are  about  two  orders  of  magnitude  below  the  concentration  of  C3-  Wachi,  in  a  later 
paper  (Reference  20)  reports  a  C3  thermochemical  model  which  exhibits  agreement  from  both  the  second 
law  and  third  law  analyses  and  also  with  the  theoretical  C3  model  of  Strauss  and  Thlole  (Reference  26). 

This  model  is  compared  with  the  JANAF  values  and  other  experimental  and  theoretical  results  in  Table  1. 

The  difference  between  the  two  C3  models  is  slight  as  shown  by  the  carbon  phase  diagram  in  Figure  1. 

These  results  certainly  provide  a  valid  experimental  basis  for  the  JANAF  carbon  vapor  thermochemistry. 

One  must  bear  in  mind  however  that  the  JANAF  predicted  triple  point  (P  ~  100  atm;  T  -  7400°R)  exhibits 
poor  agreement  with  reported  triple  point  states  (References  10,  12,  13).  This  controversy  is  re¬ 
viewed  in  the  next  section. 

2.2  CARBON  MELTING  AND  CONDENSED  PHASE  POLYMORPHIC  FORMS 

Konig  (Reference  18)  in  his  review  of  carbon  melting  states  that  "although  the  melting  tempera¬ 
tures  reported  by  various  authors,  with  few  exceptions,  are  in  good  agreement,  the  question  as  to 
whether  carbon  is  meltable  at  atmospheric  pressure  has  not  yet  been  answered  satisfactorily." 
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Thermochemical  properties  for 
Cj ,  C-j.  C<t  and  C5  Identical 

for  both  predictions 


Temperature  -®R 


Variation  In  the  carbon  phase  diagram  with  C,,  species  thermochemistry 


Konig  stated  this  observation  in  1947  and  in  the  nearly  thirty  years  hence  the  uncertainty  regard¬ 
ing  the  melting  response  of  carbon  is  not  much  improved.  Many  researchers  have  studied  the  high 
temperature  and  pressure  phase  equilibrium  of  carbon  with  the  intent  of  locating  the  triple  point. 

These  results  are  summarized  in  Table  2.  The  high  pressure  triple  point  data  in  Table  2  below  7400°R 
are  felt  to  be  of  doubtful  validity  for  the  following  reasons. 

•  Oata  below  7400nR  have  generally  been  derived  from  resistively  heated  specimen  which 
exhibit  internal  temperature  profiles  which  peak  in  the  center.  Since  the  melting  occurs 
at  the  specimen's  center,  and  the  measured  temperatures  are  surface  or  near  surface  values, 
they  are  probably  low. 

•  No  allowance  was  generally  made  for  carbon  vapor  absorption  of  the  surface  radiative  flux, 
again  causing  the  measured  temperature  to  be  below  the  actual  surface  temperature. 

Most  notable  among  the  recent  work  in  this  area  is  that  of  Schoessow  (Reference  10),  Oiaconis, 
et  al . ,  (Reference  27)  and  Whittaker,  et  al . ,  (Reference  19). 

Schoessow's  experimental  technique  reduced  the  effects  of  carbon  vapor  absorption  on  temper¬ 
ature  measurements  by  defining  a  power/temperature  correlation  for  the  black  body  hole  in  the  resis¬ 
tively  heated  specimen.  The  correlation  was  defined  at  the  lower  power/temperature  levels  where  carbon 
vapor  effects  are  minimal.  By  extrapolating  the  correlation  to  higher  power  levels,  the  melt  temper¬ 
ature  was  evaluated  from  the  input  power  at  melt.  Two  potential  effects  which  would  cause  the  measured 
temperatures  to  be  lower  than  actual  are: 

•  The  large  size  of  the  internal  black  body  cavity  and  geometry  of  the  specimen  so  as  to 
allow  substantial  thermal  gradients  to  exist  within  the  cavity.  Such  thermal  gradients 
would  tend  to  reduce  the  measured  black  body  temperature. 

•  Although  precautions  were  taken,  use  of  a  disappearing  filament  pyrometer  for  black  body 
temperature  measurements  is  such  that  any  radiation  absorption  by  the  carbon  vapor  in  the 
cavity  would  be  maximized  because  of  the  wide  wavelength  band  of  these  instruments. 

Although  it  is  difficult  to  assign  quantitative  values  to  these  uncertainties  in  Schoessow's 
temperature  data,  one  concludes  that  these  data  are  probably  the  best  available. 

The  arc  heating  technique  employed  by  Oiaconis,  et  al.,  (Reference  27)  to  measure  the  melt  tem¬ 
perature  of  ATJ-S  and  pyrolytic  graphite  proved  reasonably  successful  and  substantiates  the  data  of 
Schoessow.  The  primary  uncertainty  in  this  experimental  technique  results  from  the  uncertainty  in  the 
surface  emissivity.  By  measuring  the ’back  surface  of  the  ablating  specimen,  the  carbon  vapor  absorp¬ 
tion  effects  were  minimized,  although  uncertainties  resulting  from  temperature  gradients  within  the 
material  were  introduced. 
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Van  Vec hli'ii  (Reference  28)  has  recently  developed  a  scaling  theory  for  predicting  the  melting 

I „j  ,.„valent  crystal'.,  lie  predicts  a  i.irhon  melt  temperature  oi  «746"«  which  is  suhstan 

tially  above  any  of  the  experimental  measurements.  This  predicted  melt  temperature  is  about  600  R 
above  the  highest  measured  graphite  ablation  temperatures  in  both  the  AEOC  aeroball istic  range  (Ref¬ 
erence  29)  and  the  AFFOL  50  MW  arc  (Reference  30).  Graphite  and  carbon/carbon  ablation  temperatures 
in  the  high  pressure  regime  (Pt  -100  atm)  in  both  facilities  are  in  the  range  of  7800°R  -  8100°R. 

Whittaker's  (Reference  19)  experiments  with  laser  heated  spinning  cylindrical  samples  of  pyro¬ 
lytic  graphite  have  resulted  in  a  unique  set  of  carbon  melting  response  data.  Whittaker  observes  at 
a  total  chamber  pressure  of  ~0.2  atmospheres  and  a  temperature  of  6800°R  ±  54°R  a  liquid  phase  on  the 
surface  which  can  be  observed  flying  from  the  spinning  rod.  Whittaker,  however,  does  estimate  the 
local  vapor  pressure  to  be  150  times  the  static  pressure  which  implies  a  local  pressure  of  -30  atm. 

Smaller  droplets  thrown  off  by  the  spinning  rod  are  observed  to  maintain  their  integrity  while  vapor¬ 
izing,  whereas  larger  droplets  break-up  into  20-30  fragments  about  0.3  msec  after  leaving  the  sample. 
Whittaker  has  captured  the  carbon  matter  leaving  the  spinning  rods  on  copper  slabs.  At  static  pres¬ 
sures  above  0.2  atmospheres,  the  material  deposited  on  the  slab  has  a  uniform  black  appearance.  When 
tins  material  is  removed  from  the  surface  of  the  slab  and  the  surface  is  examined  it  exhibits  a  cratered 
appearance.  These  craters  are  felt  to  result  from  bombardment  of  the  copper  slab  with  globs  of 
liquid  carbon.  At  pressures  below  0.2  atmospheres  the  deposit  has  a  distinct  silvery  appearance  and 
a  considerably  different  character.  It  consists  of  a  relatively  dense  sheet  that  breaks  into  strips 
which  consist  mainly  of  chaoite,  a  carbyne,  and  several  other  linear  forms  of  polymorphic  carbon. 

The  linear  transitional  polymorphic  forms  of  carbon  identified  by  Whittaker  provide  heretofore 
unavailable  insight  into  the  high  temperature  response  of  the  carbon  crystalline  structure.  Whittaker 
has  observed  that  the  solid  phase  transformations  are  slow  at  low  temperature,  and  the  rate  increases 
as  the  triple  po'it  is  approached.  Whittaker  believes  the  linear  carbon  polymorphs  are  confined  to  the 
temperature  regime  of  4600°R  -  6850°R.  To  date  eight  polymorphic  forms  are  known,  and  it  is  felt  there 
may  be  others  yet  unidentified.  Evidence  exists  which  indicate  that  graphite  transforms  to  chaoite  at 
4700°R  and  that  it  is  6  carbyne  which  melts  to  produce  liquid  carbon  at  -6850°R. 

It  is  apparent  that  Whittaker's  results  are  in  conflict  with  the  results  of  Schoessow  and  others 
who  have  placed  the  triple  point  pressure  at  -100  atm,  although  there  is  some  uncertainty  regarding  the 
local  carbon  vapor  pressure  in  his  tests.  In  addition,  the  linear  polymorphic  forms  of  carbon  above 
4700°R  identified  by  Whittaker  are  a  unique  finding.  Whittaker  is  quick  to  point  out  that  these  linear 
forms  are  extremely  unstable  upon  cooldown  and  therefore  are  not  readily  available  to  be  identified  on  a 
post-test  specimen.  Cooldown  rates  of  3500°R/sec  must  be  increased  to  ~50,000°R/sec  to  freeze  the  linear 
polymorphic  forms. 
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The  review  of  work  relation  to  carbon  melting  and  the  high  temperature  ana  pressure  regime  of 
the  carbon  phase  diagram  is  presented  to  illustrate  the  uncertainty  which  currently  exists  in  estab¬ 
lishing  the  high  temperature  thennochemical  response  of  graphitic  materials.  Whittaker's  results 
which  show  an  unstable  met  formation  at  6850"M  must  be  rationalized  with  the  high  pressure  graphite 
ablation  data  which  consistently  show  surface  brightness  temperatures  in  the  range  7800°R  -  8100  R 
Possibly  the  rates  controlling  the  transformation  to  the  various  linear  polymorphic  carbon  states 
are  such  as  to  allow  the  material  to  be  heated  above  the  equilibrium  melt  temperature.  It  is  obvious 
that  questions  regarding  the  mechanisms  of  carbon  melt  will  persist  for  some  time.  At  the  present  time, 
however,  the  ablation  test  data  must  be  relied  on  for  predictions  of  carbon  and  graphite  ablation  dur¬ 
ing  reentry. 

2.3  CARBON  SUBLIMATION  RATE  UNCERTAINTIES 

This  review  of  graphite  sublimation,  the  rate  physics,  and  the  current  modeling  techniques  is 
presented  in  four  subsections.  Section  2.3.1  presents  an  overview  of  the  mass  and  energy  transfer 
models  used  to  predict  the  thermochemical  ablation  response  of  graphitic  materials.  Section  2.3.2 
reviews  background  and  development  of  carbon  kinetic  sublimation  models.  Section  2.3.3  addresses  the 
uncertainty  associated  with  the  state  of  the  carbon  vapor  species  adjacent  to  the  ablating  material 
when  sublimation  is  kinotically  controlled.  Some  models  assume  the  gaseous  carbon  vapor  species  equil¬ 
ibrate  while  other  models  assume  no  gas  phase  reactions.  Finally,  the  experimental  verification  of  sub¬ 
limation  kinetics  in  the  intermediate  temperature  regime  is  reviewed  in  Section  2.3.4. 

2.3.1  Mass  and  Energy  Transfer  Models  for  Graphite  Ablation 

This  section  presents  a  brief  review  of  the  heat  and  mass  transfer  models  and  assumptions  implicit 
in  the  thennochemical  ablation  predictions  subsequently  presented.  The  differential  equations  used  to 
describe  the  thermochemical  ablation  of  graphitic  materials  reviewed  herein  are  based  on  the  hypersonic 
approximation  of  lees  (Reference  31).  Further  assuming  boundary  layer  similarity  and  unity  Lewis-Semenov 
number,  results  in  a  direct  solution  of  the  mass  and  energy  balance  relations.  Details  of  this  reduction 
are  outl ined  below. 

Mass  Transfer 

The  mass  loss  from  a  graphitic  surface  (assuming  no  particulate  mass  loss)  is  directly  related 
to  the  surface  recession  rate  by  the  simple  expression 

ms  =  (pu)w  =  ps  (1) 
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which  is  illustrated  by  the  mass  balance  in  Figure  2. 
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Figure  2.  Thermochemical  mass  balance  for  an  ablating  graphitic  material. 

For  graphite  ablation  in  air  tne  mass  conservation  is  described  by  the  following  set  of  equations 


KN  +Ko  +KC  =1 

WWW 

(2) 

Wcw  =  P*  =  pi\  - 

(3) 

(pu)  =  0  =  psK  -  (pD^L 

w  w  0  J 

(4) 

(p n)  n  =  0  =  ,,sKn  ■  (pD^.L 

w  w  v  J 

(5) 

The  gradients  of  the  individual  chemical  elements  at  the  wall  are  given  by  an  equation  of  the 
form,  e.g. ,  carbon. 


(6) 


where  for  unity  Lewis-Semenov  number  the  convective  heat  transfer  to  the  wall  by  conduction  and  dif¬ 
fusion,  based  on  a  Stanton  number  becomes 


xonv  ■ 


■(£$).■  '.Vh  K  - H») 


(7) 


Thus,  combining  equations  (6)  and  (7)  to  solve  for  (3Ki/3y)u  in  terms  of  the  enthalpy  potential  and  the 
Stanton  number  yields  an  expression  for  any  of  the  elmental  constituents  at  the  wall 


+ 


.e* 


-jsi. _ \ 

peUeCH/ 


(8) 


where 


B' 


=  _£i 


p  u  cu 

Ke  e  H 


and  the  following  boundary  conditions  are  imposed  for  carbon  ablation  in  air, 
Kce  =  0  Kne  -  0.768  =  0.232  Kc(s)  =  1 
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F c|ua t i on  (8)  shows  that  the  mass  fraction  of  the  element  carbon  at  the  wall  is  determined  solely  by  the 
mass  flux  of  elemental  carbon  injected  at  the  surface  as  vapor.  Likewise,  equation  (8)  specifies  the 
elemental  mass  fraction  of  both  oxygen  and  nitrogen  at  the  wall  in  terms  of  the  nondimensional  mass 
transfer  parameter  B1.  A  stoichiometric  expression  then  directly  relates  the  elemental  mass  fractions 
to  the  chemical  composition  of  the  ablation  species  at  the  surface. 


u.M. 


(9) 


With  the  mass  transfer  equations  complete  the  classical  B'/T  ablation  curves  can  be  specified  by  defin 
ing  the  thermochemical  state  of  the  ablation  products.  Much  flexibility  exists  when  specifying  the 
thermochemical  state  of  the  ablation  species,  such  as 


1.  chemical  equilibrium 

2.  kinetic  controlled  sublimation  (which  imposes  an  additional  constraint  on  the  carbon  vapor 
species  surface  mass  transfer),  frozen  carbon  vapor 

3.  kinetic  controlled  sublimation,  equilibrium  carbon  vapor. 

Examples  of  B'/T  ablation  curves  are  presented  in  Figure  3,  showing  the  variation  in  the  predicted  sur¬ 
face  temperature  with  variations  in  the  assumed  thermochemistry.  To  predict  the  graphite  ablation  re¬ 
sponse  a  surface  energy  balance  must  be  coupled  to  both  the  mass  transfer  equations  and  the  in-depth 
heat  conduction.  The  energy  balance  expressions  are  developed  subsequently. 


Energy  Transfer 

The  thermochemical  energy  balance  equation  for  an  ablating  graphite  surface  is  given  below. 


p  U  Cu 
'  e  e  H 


‘  »„)  *  V.dj„ 


-  -  (pu)wHw  +  msHs  •  ^cond  =  0 


(10) 


where 


H  =  E.K.H  . 
w  J  j  J 


(n; 


This  surface  energy  balance  is  shown  in  Figure  4 
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Figure  4  Thermochemical  energy  balance  for  an  ablating  graphite  material . 
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The  blowing  to  nonblowing  heat  and  mass  transfer  ratio  for  laminar  flow  over  an  ablating  graphite 
surface  is  best  correlated  y  (Reference  32) 


CH 


2AB'o  _ 

(gZAB-o  .,) 


where  A  =  0.7 


(121 


Fquation  10,  when  coupled  to  a  standard  finite  difference  solution  to  the  indepth  conduction 
equation  enables  a  description  of  the  thermochemical  ablation  response  of  a  graphitic  surface.  Specific 
modeling  techniques  for  predicting  the  kinetic  sublimation  response  of  a  graphitic  surface  are  reviewed 
in  the  next  section. 

2.3.2  Background  and  Modeling  of  Sublimation  Kinetics 

Graphite  sublimation  is  kinetically  controlled  under  free-molecular  flow  conditions  and  vapor¬ 
ization  coefficients  have  been  measured  at  these  low  temperature  and  pressure  sublimation  states 
(Reference  7).  Before  further  consideration  is  given  to  the  sublimation  response  in  the  free  molecular 
regime,  one  must  establish  the  need  for  modeling  sublimation  kinetics  in  the  high  pressure/highly  con¬ 
vective  regime  of  interest.  To  demonstrate  how  the  existence  of  sublimation  kinetics  can  be  established, 
the  process  first  must  be  described. 

Kinetic  controlled  sublimation  is  modeled  by  using  a  Knudsen-Langmuir  expression  of  the  form. 


which  requires  that  Ihe  ith  carbon  vapor  species  partial  pressure  is  always  below  the  equilibrium  value 
for  a  non-zero  mass  flux.  Because  the  species  partial  pressures  at  the  wall  must  be  below  the  equilib¬ 
rium  values,  the  rate  controlled  sublimation  temperature  must  always  be  above  the  equilibrium  diffusion 
controlled  value.  The  fact  that  recent  high  quality  graphite  ablation  temperatures  are  above  the  values 
predicted  by  either  the  Dolton  or  Kratsch  equilibrium  diffusion  controlled  ablation  models  suggests 
that  either  the  diffusion  controlled  assumption  is  invalid  or  the  thermochemical  data  are  incorrect. 

The  early  modeling  of  carbon  sublimation  kinetics  was  done  by  Dolton,  et  al . ,  (Reference  22). 

This  initial  model  was  simplified  by  limiting  ablation  products  to  the  following  species  N2>  U2>  CO,  CN, 
C2N2’  VC16  Where  the  Production  of  CN  was  assumed  to  be  equilibrium  controlled  and  CO  production  was 
equi  1  i brium/di f fusio.i  controlled.  Carbon  sublimation  was  controlled  by  the  Knudsen-Langmuir  rate  equa¬ 
tions  shown  above.  In  addition  the  carbon  vapor  species  were  assumed  to  not  react  in  the  gas  phase 
(i.e.,  the  carbon  vapor  concentrations  were  controlled  by  the  sublimation  process).  Application  of 
this  model  illustrated  the  controlling  effects  which  sublimation  kinetics  have  on  the  predicted  graphite 
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ablation  response.  Sensitivity  studies  were  even  made  to  assess  the  significance  of  uncertainties  in 
the  species  vaporization  coefficients  ( t tc . )  on  the  ablation  predictions. 


Siih'.ei|iii'iit  ipapliite  kinetic  sublimation  models  developed  by  /orliiq  (Holer;1 'nee  33)  and  Kralsdi, 
id  al.,  (Reference  14)  were1  less  restrictive  with  respect  to  the  number  of  ablation  products  (CN,  L,^, 
Cj,N)  and  the  reactivity  of  the  carbon  vapor  species  (i.e.,  equilibrium  carbon  vapor).  Both  models  are 
considered  an  improvement  over  the  initial  model i ny  of  Dolton,  et  al.  Recent  application  of  the 
Aerotherm  Chanical  Equilibrium  (ACE)  computer  code  (Reference  35)  to  graphite  kinetic  sublimation  pre¬ 
dictions  has  further  advanced  the  generality  of  these  calculations.  No  restriction  exists  with  respect 
to  the  number  of  ablation  products  and  the  sublimation  products  can  either  be  assumed  to  equilibrate 
or  be  frozen.  Application  of  this  ablation  model  has  provided  some  insight  into  the  significance  of 
the  assumed  carbon  vapor  state  at  the  surface  (Reference  36).  These  results  are  reviewed  in  the  sub¬ 
sequent  section. 

2.3.3  Effect  of  th e  Surface  Carbon  Vapor  State 

Graphite  ablation  predictions  in  Reference  36  illustrate  the  sensitivity  of  kinetic  sublimation 
predictions  to  the  assumed  state  of  the  carbon  vapor  species  at  the  surface.  Results  of  these  calcula¬ 
tions  show  that  at  pressures  below  one  atmosphere  the  shift  to  higher  temperatures  with  increasing 
f’eUeCm  ^at  constant  B<)  relatively  insensitive  to  the  state  of  the  carbon  vapor  species  at  the  sur¬ 
face.  This  result  is  shown  in  Figure  5.  However,  at  pressures  above  10  atmospheres,  the  sensitivity 
of  B' /temperature  curves  to  the  assumed  carbon  vapor  state  at  the  surface  increases  substantially . 

This  result  is  shown  in  Figure  5  at  100  atmospheres.  In  other  words,  at  high  pressures,  if  carbon 
sublimation  is  kinetically  controlled,  but  the  vapor  species  equilibrate  in  the  gas  phase  the  theo¬ 
retical  B '/ temperature  curve  approaches  the  equilibrium  diffusion  controlled  curve,  as  shown  in  Figure 
5.  Thus,  the  low  to  moderate  pressure  sublimation  tests  in  the  NASA  Ames  AEHS  facility  will  not  9e 
sensitive  to  uncertainties  in  the  carbon  vapor  state  at  the  ablating  surface.  A  more  complete  dis¬ 
cussion  of  the  carbon  vapor  state  uncertainties  and  implications  thereof  as  this  relates  to  ablation 
predictions  for  flight  is  presented  in  Section  6. 

2.3.4  Experimental  Identification  of  Vaporization  Kinetics  Effects 

It  is  apparent  from  the  discussions  in  the  previous  two  sections  that  kinetically  controlled 
sublimation  has  the  effect  of  reducing  the  sublimation  rate  from  an  ablating  surface  under  a  fixed 
set  of  environmental  conditions.  This  result  will  cause:  (1)  the  surface  ablation  temperature  to  rise 
and  (2)  the  surface  recession  rate  to  decrease.  These  effects  depend  on  the  convective  mass  transfer  rate 
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within  the  boundary  layer.  This  behavior  is  evident  in  the  comparison  of  equilibrium  diffusion  con¬ 
trolled  and  kinetically  controlled  B1  curves  shown  in  Figures  3  and  5.  This  section  is  designed  to 
show  how  the  kinetic  sublimation  shift  to  higher  temperatures  varies  with  the  convective  hfat  and 
mass  transfer  rate  within  the  boundary  layer.  Physically  one  would  expect  the  kinetic  sublimation 
predictions  to  agree  with  the  equilibrium-diffusion  controlled  predictions  when  the  mass  transfer  rate 
approaches  zero.  The  analysis  developed  below  illustrates  this  result  by  considering  the  sublimation 
kinetics  and  boundary  layer  diffusion  to  be  two  resistances  in  series.  Such  an  approach  yields  the 
following  implicit  expression  for  the  ratio  B1 rate^'equi  1  ibri urn  (Reference 
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where 


a  ~  vaporization  coefficient 
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Since  11 1  ....  .  z  fn  (T,P)  this  expression  shows  that  for  the  two  limiting  conditions  of 

'  equil  ibrium 


U  cM->  0  and  i,  U  CM-»  -  the  following  values  of  B1  rate/B 'equii  ibrium  are  attained' 


eue''M~'  u  “”'J  ‘e  e  M  - ”  “  rate'  equi 

B'rate 


Limit 

peUeCM  0  B  equilibrium 


=  1 


B1 


rate 


Limit 

peUeCM  *  °°  B  equilibrium 


=  0 


These  limits  are  readily  derived  from  expanded  forms  of  the  above  expressions. 

B1 


Limit 

p  U  C. 
ee  e  1 


/B  rate  \  peUeCM 

’M  *  °  \^J  W 


eUeCM  „  .  B  rate  .  , 
K  +  B7"  "  ‘  1 


equil 


F 


rate  .  1 


equil 


Limit 

p  U  C 
1  e  e 


B  rate  _  1  j . 

J(  ap  V  + 

4au 

-*■  00 

*1^17  2| 

VV7SV 

p  U  CUK 
e  e  M 

p  U  CUK 
e  e  M 

rate  _  „ 
p - 0 


equil 


17 


TeM?e.eATO£G  —  e- 

(a)  1.0  Atmosphere 

Figure  5.  Sensitivity  of  kinetic  sublimation  predictions  to  the  carbon  vaoor  state 


(b)  TOO  atmospheres 


Thus,  properly  designed  experiments  which  cover  a  large  range  in  the  value  of  the  convective  mass 
transfer  coefficient  n  U  C,.  can  establish  the  existence  of  rate  controlled  sublimation  response  for 
ablation  conditions  of  interest.  Figure  3  shows  that  if  the  kinetic  sublimation  model  recommended 
by  Dolton,  et  al.,  (Reference  22)  using  the  JANAF  data  is  assumed,  a  factor  of  5  variation  in 


'e^e^M  's  reRu'ret*  t0  Produce  a  200"R  shift  in  temperature  at  0.1  atmosphere  for  B'  1.5.  To 
achieve  a  factor  of  5  variation  in  p  U  CM  in  the  staunation  region  of  a  blunt  cyclindrical  model  by 


geometric  variations,  the  model  diameter  would  have  to  vary  by  about  a  factor  of  25.  It  is  immedi¬ 


ately  obvious  that  this  is  not  possible  within  the  constraints  of  most  test  facilities.  Therefore, 


the  best  means  of  achieving  the  desired  variation  in  is  by  varying  the  velocity  of  the  con¬ 


vective  flow  from  subsonic  to  supersonic. 


In  order  to  maintain  the  surface  temperature  at  sufficiently  high  values  to  sustain  subli¬ 
mation  while  varying  the  convection  energy  and  mass  transfer  coefficients,  an  additional  source 
of  energy  flux  into  the  surface  must  be  available.  Two  means  of  achieving  this  are: 

•  Incident  radiation  combined  with  convective  heating 
a  Resistively  heated  sample  with  convective  heating 

Both  means  of  conducting  such  ablation  experiments  were  thoroughly  studied.  It  was  established  that 
the  combined  radiative/convective  technique  is  superior  to  the  combined  resistive/convective  technique. 

Properly  designed  radiativo/convective  heating  tests  can  provide  a  sufficiently  large  variation 
in  the  convective  coefficient  while  maintaining  the  ablation  state  at  sufficiently  high  B'  values  in 
the  sublimation  regime.  In  fact,  the  equilibrium  diffusion  limited  condition  ((>eUeCM  ~0)  is  approached 
at  one  atmosphere  test  conditions  with  only  radiative  heating.  The  results  of  Davy  and  Bar-Nun 
(Reference  38)  substantiate  this  as  a  viable  technique  for  measuring  the  equilibrium  vaporization 
temperature  at  pressures  above  0.1  atmosphere,  although  free  convection  currents  and  carbon  vapor 
condensation  tend  to  complicate  the  interpretation  of  the  results  at  the  one  atmosphere  test  condition. 


The  only  combined  radiative/convective  facility  currently  in  existence  with  the  required  flex¬ 
ibility  for  such  a  test  program  is  the  Advanced  Entry  Heating  Simulator  (AEHS)  at  Ames  Research  Center 
(Reference  39).  This  facility  was  used  by  Davy  and  Bar-Nun  to  conduct  "radiative-only"  ATJ-S  ablation 
experiments  and  by  Wakefield  and  Peterson  for  combined  radiative/convective  ATJ-S  ablation  tests.  The 
ablation  test  program  described  in  Section  3  has  extended  the  work  of  these  investigators  by  combining 
high  enthalpy  supersonic  and  subsonic  flow  conditions  with  the  radiative  heat  flux  maintaining  the 
ablation  response  at  sufficiently  high  B'  values  to  observe  the  potential  temperature  shifts  with  ,->e Ue CM . 
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2.4  PARTICULATE  MASS  LOSS  UNCERTAINTIES 

Accurate  interpretation  of  ablation  test  data  depends  on  one's  ability  to  accurately  assess 
the  state  of  matter  relieved  during  the  ablation  process.  This  question  has  been  addressed  by  nearly 
all  researchers  who  have  worked  with  the  hyperthermal  response  of  graphitic  materials.  The  Question 
also  relates  closely  to  the  uncertainty  of  the  high  temperature  phase  equilibrium  of  carbon. 

Lundell  and  Dickey  (Reference  24)  devised  a  technique  to  photograph  the  particulate  mass  loss 
from  ATJ  graphite  models  tested  in  the  NASA  Ames  Heat  Transfer  Tunnel.  They  estimate  from  their 
results  that  the  particulate  flux  becomes  a  significant  component  of  the  total  mass  loss  at  tempera¬ 
tures  above  6650°R.  Because  of  the  two  phase  nature  of  ATJ  graphite  (i.e.,  filler  particles  held 
together  in  a  binder  matrix)  it  was  rationalized  that  preferential  ablation  of  the  binder  material 
allowed  filler  particles  to  be  released  and  thus  removed  in  the  solid  state. 

Lundell  and  Dickey  performed  a  similar  set  of  ablation  experiments  using  over  80  different 
types  of  graphitic  materials  with  the  intent  of  evaluating  the  role  of  the  material  microstructure 
in  controlling  particulate  mass  loss  (Reference  25).  This  second  series  of  ablation  tests  was  sim¬ 
ilar  to  the  previous  series  in  all  respects.  The  graphitic  materials  tested  included  commercial  and 
developmental  grades  of  artificial  graphites,  both  two  and  three  dimensional  carbon/carbon  composites 
seeded  with  refractory  compounds,  and  several  special  materials  such  as  pyrolytic  graphite,  measophase 
graphite,  glass  carbon,  and  natural  graphite.  Results  of  these  tests  showed 

1.  The  extreme  variation  in  these  materials'  microstructure  has  little  effect  on  their  ablation 
performance  at  these  low  pressure  (-4  atm)  test  conditions. 

2.  The  observed  particulate  mass  loss  from  single  phase  carbon  materials  indicates  preferen¬ 
tial  ablation  is  not  the  only  explanation  for  particulate  mass  loss. 

3.  No  correlation  relating  the  microstructure  of  commercially  available  graphites  to  the  par¬ 
ticulate  mass  loss  was  evident  from  these  tests. 

These  results  of  Lundell  and  Dickey  show  that  for  the  low  to  moderate  pressure  test  conditions  consid¬ 
ered,  particulate  mass  loss  from  graphitic  materials  does  exist,  but  does  not  correlate  with  the  ma¬ 
terials'  microstructure.  One  must  keep  in  mind,  however  that  the  particulate  mass  loss  observed  by 
lundell  may  be  unique  of  the  test  technique  employed 

Whittaker  interprets  these  results  to  show  the  effect  of  carbon  melting  at  6850°R,  since  the 
sharp  rise  in  pa-ticulate  mass  loss  observed  by  Lundell  and  Dickey  occurs  at  a  surface  temperature  of 
6650°R.  Whittaker  and  Kintner  (Reference  40)  also  observed  particulate  mass  loss  at  temperatures 
above  6650°R  and  cite  similar  observations  by  Maahs  and  Schryer  (Reference  41)  and  Wachi  and  Gilmartin 
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(Reference  15).  It  is  only  natural  that  Whittaker  explains  the  phenomena  of  particulate  mass  loss 
in  this  temperature  regime  to  result  from  carbon  melting,  since  the  results  of  his  research  show 
a  carbon  melt  phase  to  exist. 

A  substantial  amount  of  work  was  done  under  the  Air  Force  sponsored  Passive  Nosetip  Technology 
(PANT)  Program  to  evaluate  the  effects  and  probanility  of  particulate  mass  loss  from  ATJ-S  graphite 
under  high  pressure  hyperthermal  environments  (Reference  42).  Ablation  test  data  from  both  the  AFFDL 
50  MW  arc  and  the  AEDC  Aeroballistic  Range  were  thoroughly  analyzed  based  on  the  JANAF  thermochemical 
ablation  model  with  particulate  mass  loss  and  rough  wall  convective  heating  sensitivity  studies. 

These  ablation  data  were  the  first  high  pressure  graphite  ablation  data  with  both  high  quality  surface 
temperature  and  surface  recession  data.  The  surface  temperature  and  recession  data  made  it  possible 
to  assess  the  significance  of  particulate  mass  loss  under  these  high  pressure  ablation  test  conditions. 
These  analyses  consistently  showed  particulate  mass  loss  to  be  small. 

Kratsch,  et  al.,  within  a  subtask  of  the  ASML  Exploratory  Development  Program  (EDP)  titled,  "The 
Erosion  Mechanisms  and  Improvement  of  Graphitic  Materials,"  (Reference  34),  studied  at  the  macro  and  micro- 
structural  response  of  graphitic  materials  under  hyperthermal  ablation  conditions.  A  segment  of  this 
program  specifically  studied  the  ablation  response  of  bulk  graphitic  materials  including  ATJ-S  graphite. 

This  micromechanical  erosion  model  was  applied  and  evaluated  on  the  PANT  program.  The  MDAC/EDP 
micromechanical  ablation  model  is  summarized  below. 

•  The  surface  thermochemical  ablation  response  is  kinetically  controlled. 

•  Subsurface  sublimation  occurs  by  diffusion  of  equilibrium  carbon  vapor  through  the  porous 
microstructure  in  the  near  surface  material. 

•  Subsurface  sublimation  weakens  the  surface  structure.  If  sufficiently  weakened,  the  surface 
material  can  be  mechanically  removed. 

Application  of  the  micromechanical  ablation  model  to  both  50  MW  and  ballistic  range  ablation 
data  resulted  in  no  predicted  particulate  mass  loss.  Kratsch,  et  al.,  report  that  their  use  of  the 
model  showed  particulate  mass  loss  to  generally  be  less  than  10  percent  of  the  total  with  no  particulate 
mass  loss  often  being  predicted. 

Review  of  this  work  is  included  in  this  discussion  to  illustrate  the  uncertainty  which  currently 
exists  regarding  the  existence  and  significance  of  particulate  mass  loss  from  hyperthermal  ablating 
graphite  materials.  To  quantitatively  assess  the  significance  of  particulate  mass  loss  it  must  be 
directly  measured.  To  date  no  technique  has  been  devised  to  accomplish  this  task,  thus  the  precise 
amount  of  particulate  mass  loss  cannot  be  measured.  Because  of  this  limitation  one  must  infer  the 
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SECTION  3 


EXPERIMENTAL  APPROACH 

The  experimental  approach  utilized  in  the  carbon  sublimation  study  is  described  in  this  sec- 

4 

tinn.  Section  3.1  discusses  the  test  procedure  and  approach  as  related  to  the  objectives  of  the 
study.  Section  3.?  is  a  description  of  the  tost  equipment  employed  for  the  experimental  portion  of 
the  study.  Section  3.3  reviews  the  test  program  including  both  calibration  and  sublimation  tests 
while  the  test  results  are  presented  in  Section  3.4. 

3. 1  TEST  PROCEDURE/APPROACH 

The  review  of  the  Kriudsen-Lamjmuir  kinetic  sublimation  model  in  Section  2.3  shows  that  sub¬ 
limation  reactions  which  are  ki net i ca 1 1 y  controlled  exhibit  a  dependence  on  the  convective  mass 
transfer  coefficient  (i'euec|v|)-  Because  of  this  dependence  the  best  way  to  assess  the  existence  and 
significance  of  sublimation  kinetics  is  to  perform  carbon  sublimation  tests  under  varying  convective 
conditions.  In  order  to  maintain  ablation  in  the  Sublimation  regime,  some  means  of  augmenting  the 
convective  heat  input  (such  as  incident  radiation  or  Inductive  heating)  is  required  as  discussed  in 
Section  2.3 

A  review  of  available  test  facilities  with  multiple  heating  modes  identified  the  NASA  Ames 
Advanced  Entry  Heating  Simulator  (AEHS)  as  the  optimum  facility  for  this  test  series.  The  AEHS 
facility  has  the  capability  to  subject  a  test  specimen  to  simultaneous  convective  and  radiative 
heating.  The  peak  radiation  flux  which  can  be  generated  with  its  argon  driven  arc  imaging  radiation 
source  is  approximately  2500  Btu/ft?  sec  which  is  adequate  to  maintain  carbon  in  its  sublimation 
state  under  a  near  zero  convective  heating  environment. 

Therefore  the  ablation  test  procedure  used  on  this  program  was  to  expose  models  to  two  nomi¬ 
nal  pressure  conditions  (i.e.,  0.1  and  -  1.0  atmospheres)  and  various  convective  conditions. 

The  variation  in  the  convective  heat  and  mass  transfer  rate  was  achieved  by  varying  the  arc  heated 
air  velocity  over  the  models.  This  procedure  was  used  since  it  was  impossible  to  vary  the  model 
size  sufficiently  to  gain  the  desired  variation  in  PeueCfl).  The  size  of  the  radiation  beam  with  a 
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90  percent  uniformity  is  shown  in  Figure  6  to  be  0.5  inch  in  diameter  at  a  0.25  inch  standoff. 

Thus,  the  0.5  inch  diameter  corresponds  to  an  upper  limit  on  the  model  size.  Since  the  pyrometer 
which  views  the  ablating  surface  is  located  -  30  inches  from  the  model,  the  model  diameter  can't  be 
much  smaller  than  0.25  inch  with  any  assurance  of  measuring  accurate  surface  temperature  data. 

Thus,  the  practical  size  range  available  for  models  lies  between  0.25  and  0.5  inch.  A  factor  of 
two  variation  in  nose  radius  corresponds  to  a  40  percent  variation  in  peugCM  which  is  insufficient 
for  studying  sublimation  kinetics  as  discussed  in  Section  2.3.4. 

The  flow  over  the  models  in  this  test  series  was  varied  from  low  subsonic  velocities  (-  100 
ft/sec)  to  supersonic  flow.  The  variation  in  PgUgC^  at  the  ~  0.1  atmosphere  and  -  1.0  atmosphere 
test  conditions  was  approximately  a  factor  of  ~  40  and  ~  30  respectively.  Since  the  supersonic 
flow  condition  controls  the  maximum  temperature  shift  due  to  kinetic  sublimation  response  with  a 
fixed  configuration,  a  0.5-inch  diameter  model  at  0.1  and  1.0  atmosphere  test  conditions  experiences 
peak  hear  and  mass  transfer  rates  of  0.09  and  ~  0.19  lbm/ft2sec  respectively.  Thus  the  predictions 
in  Figure  3  show  that  at  a  nominal  B1  value  of  unity,  the  maximum  shift  in  surface  temperature  antic¬ 
ipated  at  the  •  0.1  and  -  1.0  atmosphere  conditions  are  -  220°R  and  -  120°R,  respectively.  It  was 
known  from  the  inception  of  this  program  that  these  potential  temperature  shifts  were  close  to  the 
anticipated  >2  percent  uncertainty  in  the  pyrometer  daca. 

3.2  TEST  EQUIPMENT 

This  section  describes  the  facility  hardware,  equipment  modifications  required  for  this  pro¬ 
gram,  and  the  instrumentation  used  in  data  acquisition.  The  ablation  model  configuration  is  also 
presented  and  the  procedures  used  to  establish  the  specific  model  design  employed  for  this  test 
series  are  discussed.  The  AEHS  facility  is  described  in  Section  3.2.1,  the  facility  modifications 
are  discussed  in  Section  3.2.2,  the  instrumentation  and  models  are  described  in  Sections  3.2.3  and 
3.2.4,  respectively. 

3.2.1  AEHS  Facility 

The  experimental  portion  ol  this  program  was  performed  in  the  Advanced  Entry  Heating  Simu¬ 
lator  (AEHS)  facility  at  the  NASA  Ames  Research  Center.  The  AFHS  facility  is  shown  in  Figure  7 
and  described  in  more  detail  in  Reference  39.  It  consists  of  four  main  subsystems:  the  convective 
source,  the  radiant  source,  test  chamber/vacuum  system  and  the  test  model  support  and  insertion 
system. 

The  convective  heating  source  is  provided  by  an  arc  heater  coaxially  mounted  on  the  test 
chamber/radiant  source  centerline  (see  Figure  7).  Several  arc  heaters  are  available  for  use  in  this 
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Figure  7.  Advanced  entry  heating  simulator. 


facility  and  consist  of  both  the  Linde  type  (i.e.,  N-4001)  and  the  Ames  developed  constricted  type. 
The  typical  convective  test  stream  used  in  this  facility  is  low  pressure,  high  enthalpy,  supersonic 
flow  providing  stagnation  convective  heating  rates  of  more  than  2500  Btu/ft2  sec  with  a  nominal 
model  size  of  0.5  inch  diameter.  In  order  to  meet  the  objectives  of  the  carbon  sublimation  study, 
it  was  necessary  to  modify  the  constricted  arc  heaters.  These  modifications  are  discussed  in  the 
following  section.  The  arc  heaters  are  electrically  powered  by  a  3  MW  dc  power  supply  consisting 
of  8  rectifier  modules,  each  capable  of  600  amp  output  with  an  open  circuit  voltage  of  1125  volts. 
Simulation  of  planetary  atmosphere  (air,  N2,  He,  C02,  AR)  is  provided  by  a  multiple  gas  mixing  and 
metering  system  capable  of  simultaneously  injecting  three  different  gases  over  a  100:1  mass-flow 
range  of  1.1  1  bin/ sec  to  0.011  lbm/sec.  The  dc  arc  current  and  mass-flow  rates  can  be  continuously 
varied  and  maintained  during  operation.  The  following  section  describes  the  changes  in  the  meter¬ 
ing  system  necessary  for  the  lower  mass  flow  requirements  of  this  test  series. 
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The  independently  controlled  radiant  source  can  provide  incident  radiant  energy  to  the  model 
surface  while  maintaining  reasonable  target  uniformity  in  a  low  pressure  convective  environment. 

The  system  consists  of  an  Ames  developed  125  kw,  argon  vortex-stabilized  radiation  source  and 
quasi-bifocal  arc  imaging  optical  system.  The  optical  components  (source,  collector,  transfer  lens, 
and  imaging  mirror)  are  coaxial  with  respect  to  the  tunnel  centerline.  The  collector  sr.rror  cap¬ 
tures  the  output  of  the  source  transferred  arc  and  directs  the  energy  towards  the  1  meter  diameter 
imaging  mirror  via  the  transfer  lens  (see  Figure  7).  The  imaging  mirror,  constructed  around  the 
arc  heater  nozzle,  in  turn  focuses  the  radiation  onto  the  model  surface.  The  surface  contours  of 
the  mirrors  are  distorted  from  a  simple  elliptical  shape  to  affect  increased  arc  utilization  and 
radiation  distribution  uniformity.  Complete  interruption  and/or  exposure  of  radiation  is  con¬ 
trolled  manually  or  automatically  by  a  pneumatically  operated  douser  (shutter). 

The  radiant  source  spectral  characteristics  are  shown  in  Figure  8.  This  figure  shows  that 
at  0.9  |i  (the  peak  sensitivity  of  the  Thermogage  optical  pyrometer)  a  significant  amount  of  re¬ 
flected  radiation  might  be  picked  up  by  the  Thermogage  pyrometer.  To  minimize  incident  radiation 
reflectance  errors  in  pyrometer  measurements  a  chopping  device  is  used  which  periodically  occludes 
the  incident  radiative  enerqv  for  short  intervals.  The  chopper  consists  of  two  1-meter  radius 
counter-rotating  blades.  The  blades  rotate  at  20  cps  and  entirely  block  the  beam  for  3  msec  inter¬ 
vals,  during  which  time,  the  temperature  measurements  are  taken. 

The  3-meter  diameter  test  chamber  provides  a  low  pressure  nozzle  discharge  environment  and 
houses  both  the  convective  and  radiative  systems  as  well  as  the  model  support  system.  The  chamber 
is  connected  to  a  series  of  vacuum  spheres  with  a  20.129  ft3  capacity. 


A  multiple  model  support  system  provides  for  insertion  of  test  models  or  calibration  probes 
into  the  combined  convective- radio live  environment.  A  total  of  nine  models  or  probes  can  be  auto¬ 
matically  inserted,  in  any  sequence,  for  preselected  exposure  times.  Total  insertion  time  is  less 
than  1  second  including  transit  time  through  the  arc  heated,  expanded  free  jet  of  approximately 
200  msec.  The  model  struts  are  aligned  to  the  optical  axis  and  convective  stream  centerline.  A 
laser  monitored  feedback  servo-control  system  maintains  the  model  stagnation  point  at  a  fixed  loca¬ 
tion  from  the  nozzle  exit  and  within  the  focal  point  of  the  radiation  imaging  mirror. 

3.2.2  F acjj ity  Modifications 

The  basic  objective  of  the  experiments  was  to  demonstrate  the  existence  or  nonexistence  of 
carbon  Sublimation  kinetics  by  providing  large  variations  in  mass  transfer  coefficient  (Peue^) 
while  maintaining  the  surface  temperature  in  the  sublimation  regime.  The  combined  radiative/con¬ 
vective  environment  in  the  AEHS  facility  was  capable  of  maintaining  a  relatively  constant  surface 
temperature  with  variation  in  the  mass  transfer  coefficient  by  varying  the  flowrate  over  the  model. 
Pretest  analyses  showed  that  this  objective  could  be  met  by  providing  convective  flows  from  sub¬ 
sonic  tn  supersonic  with  a  gas  enthalpy  in  the  range  of  20,000  Btu/lbm  as  discussed  in  Section  3.1. 

The  arc  heater  generally  used  in  the  facility  is  the  Linde  N-4001 .  This  heater  provides 
total  enthalpy  gas  streams  of  2,000  to  20,000  Btu/lbm  (supersonic  only)  at  model  stagnation  pres¬ 
sures  to  1  atm.  The  constricted  heater,  however,  is  capable  of  gas  enthalpies  as  high  as  40,000 
Btu/lbm  at  stagnation  pressures  to  0.3  atm.  The  required  range  of  arc  heater  operation  (subsonic- 
supersonic)  combined  with  the  enthalpy  level  desired,  dictated  the  use  of  the  constrictor  heater. 
Further,  it  was  felt  certain  modifications  were  required  before  this  heater  could  run  in  a  stable 
mode  with  the  low  flow  rates  required  for  the  subsonic  cases. 

The  principal  modification  to  the  constricted  heater  was  the  replacement  of  the  existing  pin 
electrode/ballast  resistor  anode  configuration  with  a  combined  anode/nozzle  design  which  enhanced 
the  enthalpy  capability  by  providing  a  "peaked"  enthalpy  profile.  The  peaking  was  considered  neces¬ 
sary  to  provide  the  enthalpy  level  required  without  exceeding  safe  current  limits  at  the  flowrates 
required  for  the  supersonic  operating  points.  This  operating  procedure  had  the  added  benefit  of 
maintaining  a  very  low  level  of  contaminants  in  the  test  stream.  These  modifications  also  provided 
the  necessary  facility/arc  heater  interface  for  proper  mounting  to  the  test  chamber. 

The  six  arc  heater  configurations  required  for  this  test  series  are  summarized  in  Table  3. 
Configurations  3  and  6  are  the  supersonic  configurations  at  0.1  and  1.0  atm  stagnation  pressures, 
respectively.  The  enthalpy  levels  noted  for  these  configurations  assume  a  peaking  ratio  of  2.  The 


existing  facility  2.75-inch  nozzle  extension  was  augmented  with  a  new,  1.0-inch  nozzle,  thus  sig¬ 
nificantly  lowering  the  required  flowrate  to  achieve  the  1.0  atm  stagnation  pressure.  This  configu¬ 
ration  also  determined  the  longest  constrictor  required,  46  segments.  The  reduced  flowrate  for  the 
0.1  atm  stagnation  pressure  allowed  the  constrictor  to  be  shortened  to  23  segments.  The  flowrates 
noted  in  this  table  were  calculated  from  best  estimates  of  arc  heater  performance.  During  actual 
checkout/calibration  runs,  the  flowrates  were  modified  as  required  to  produce  optimum  arc  heater 
operation. 

Configurations  4  and  5  provided  the  1.0  atm,  subsonic  conditions.  Low  and  high  velocity 
streams  were  achieved  by  using  the  two  nozzles.  The  low  flowrate  noted  was  necessary  to  prevent 
sonic  conditions  at  the  nozzle  throat.  The  0.1  atm  subsonic  conditions  were  provided  by  Configu¬ 
rations  1  and  2.  For  these  conditions,  the  arc  heater  constrictor  was  shortened  to  5  segments, 
reducing  the  arc  voltage  so  that  the  arc  current  levels  could  be  maintained  at  reasonable  levels. 

In  addition,  this  increased  the  ability  of  the  very  low  flow  to  stabilize  the  arc  column,  providing 
stable  arc  operation.  As  noted  in  the  table,  the  gas  injection  rings  of  the  arc  heater  were  modi¬ 
fied  to  maintain  sonic  gas  injection  to  the  arc  chamber  over  the  large  flow  range  required. 

The  reconfigured  constricted  arc  heater  (46  segments)  is  shown  in  Figure  9.  The  design  of 
the  required  arc  heater  parts  was  done  by  Aerotherm  with  all  fabrication  provided  by  NASA  Ames.  The 
heater  parts  required  for  the  modifications  are  listed  below.  Orawings  of  these  parts  are  included 
in  Appendix  A. 


drawing  No 

Description 

7111-033 

Oownstream  Housing,  Anode 

7111-004 

Upstream  Housing,  Anode 

7111-005 

Split  Core,  0.50-2.0  Liner 

7111-006 

0.50-2.0  Liner,  Nozzle 

7111-007 

Split  Ring 

7111-009 

Split  Core,  0.5-0.75  Liner 

7111-010 

0  50-0.75  Liner,  Nozzle 

7111-011 

0  75-1.0  Nozzle  Extension 

7111-012 

Insulator,  Bolt 

7111-013 

Insulator,  Stud 

7111-014 

Outer  Insulator 

The  range  and  low  level  of  gas  flow  rates  required  for  these  arc  heater  configurations  re¬ 
quired  the  addition  of  a  gas  metering  board  to  augment  the  facility  gas  system.  The  metering  board 


Figure  9. 

Modified  NASA  Ames  constrictor  arc  heater  used  for  graphite 
sublimation  Kinetics  study  (46  segment  configuration). 


provided  the  necessary  gas  regulation  and  monitoring  equipment  to  meter  nitrogen  and  oxygen  sepa¬ 
rately,  as  required  by  tiie  constricted  arc  heater.  In  addition  necessary  valving  and  control  equip¬ 
ment  was  included  to  provide  the  argon  start  gas  to  the  heater.  The  metering  board  schematic  is 
shown  in  Figure  10.  The  design,  assembly  and  checkout  of  the  metering  board  was  accomplished  by 
Aerotherm.  The  board  was  delivered  to  Ames  and  is  now  a  permanent  part  of  the  AEHS  facility. 

3.2.3  Instrumentation 

The  following  data  were  recorded  for  the  various  test  systems: 

Arc  heater  -  voltage,  current  and  arc  chamber  pressure 

Radiation  source  —  voltage  and  current 

Convective  stream  -  heating  rate,  stagnation  pressure  and  test  cabin  pressure  history 
Model  -  surface  recession,  surface  temperature 

The  normal  procedure  for  determining  the  arc  heater  efficiency  in  facilities  of  this  type  is 
to  measure  the  energy  loss  to  the  arc  heater  wall  by  monitoring  the  cooling  water  temperature  change 
and  mass  flowrate.  With  the  heater  efficiency  known,  the  stream  bulk  enthalpy  can  be  determined 
from  the  power  input  and  gas  mass  flowrate.  The  AEHS  facility  did  not  have  this  capability  and  the 
stream  enthalpy  was,  therefore,  backed  out  from  the  measured  stagnation  point  heat  flux  through  use 
of  the  Fay  and  Riddell  stagnation  point  heat  transfer  correlation  for  the  supersonic  test  cases  and 
by  an  iteration  technique  described  in  Section  3. 4. 2.1  in  the  case  of  the  subsonic  test  points. 

All  the  above  mentioned  pressure,  heating  rate,  voltage  and  current  data,  and  pyrometer  out¬ 
put  were  recorded  on  a  high  speed  CEC  oscillograph  system  for  later  reduction  by  Aerotherm  person¬ 
nel.  The  arc  voltage,  current,  chamber  pressure,  and  test  cabin  pressure  were  also  monitored  by 
facility  gauges  and  recorded  photographically  and/or  manually  for  comparison  to  the  CEC  data.  Model 
recession  data  were  recorded  by  a  high  speed  Photo-sonics  16mm  motion  picture  camera  viewing  the 
model  in  the  horizontal  plane  183°  from  the  pyrometer  viewing  location. 

Pressure  data  were  measured  by  Statham  absolute  and  differential  pressure  transducers  of 
various  ranges  to  cover  the  different  operating  regimes  experienced  in  this  test  series. 

Heating  rate  data  were  measured  by  steady-state  circular  foil  heat  flux  gages  as  well  as  an 
Ames  manufactured  copper  slug,  capacity  type  calorimeter.  Two  types  of  steady-state  calorimeters 
were  available.  The  Hycal  was  a  0-5000  Btu/ft2  sec  calorimeter  intended  for  the  supersonic  high 
heat  flux  portion  of  the  test  series.  A  number  of  Medtherm  calorimeters  were  available  in  ranges 
of  0-1000  Btu/ft2  sec  and  0-3000  Btu/ft2  sec  for  the  lower  heating  rates.  However,  these  devices 


showed  little  repeatability  and  erratic  response  during  calibration  runs.  The  data  from  these 
dev  ices  were  ignored  and  the  llytal  and  slug  calorimeter  became  the  primary  source  of  heating  rate 

da  ta . 

Accurate  radiant  heating  rate  data  was  obtained  by  coating  the  transient  calorimeter  with  a 
high  emissivity  coating  (camphor  soot)  before  each  run.  Comparison  of  the  coated  slug  data  with 
the  steady-state  calorimeter  measurements  made  in  the  same  environment  during  the  same  run  provide 
an  in-situ  calibration  of  the  steady-state  calorimeter.  This  comparison  was  most  useful  since  the 
radiant  heating  levels  were  capable  of  altering  the  transient  calorimeter  coating  characteristics 
(i.e.  ,  removal  of  some  coating)  after  only  one  insertion.  The  steady-state  calorimeter  then  became 
the  primary  measurement  device  for  subsequent  radiant  data  during  the  same  test  run. 

The  surface  temperature  data  were  recorded  by  a  Thennogage  (S/N  422)  optical  pyrometer. 

The  instrument  viewed  the  model  through  a  barium-fluoride  window  and  an  aperture  in  the  large  "ellip¬ 
tical"  imaging  mirror.  The  viewing  angle  was  45°  off  of  the  model  axis  with  an  approximate  spot 
size  of  1/16-inch  diameter  at  a  focal  length  of  30  inches.  The  sensing  wavelength  of  this  pyrom¬ 
eter  is  centered  at  0.8  micron  and  falls  almost  entirely  within  the  band  from  0.4  to  1.1  microns. 

The  pyrometer  spectral  response  is  shown  in  Figure  11.  A  neutral  density  filter  was  used  to  extend 
the  range  of  the  pyrometer  to  the  7600°R  level. 

The  relative  size  of  the  sensing  area  and  the  model  size  required  careful  alignment  of  the 
pyrometer  to  insure  accurate  temperature  data.  To  this  end,  the  Thermogage  standard  "1 lght-source- 
replaceiient-of- the- sensor"  technique  was  used  prior  to  each  test  run.  The  location  of  the  resulting 
pyrometer  sensing  area  is  shown  in  Figure  12. 

The  model  surface  reflective  component  from  the  radiation  source  was  eliminated  from  the 
pyrometer  data  by  the  source  chopper  as  discussed  earlier  (Section  3.2.1).  The  relatively  fast 
response  time  of  the  Thermogage  pyrometer  and  associated  electronics  (on  the  order  of  ~  3  msec) 
allow  the  surface  temperature  data  to  be  obtained  during  this  "blocked"  interval  during  each  chopper 
cycle.  An  example  of  a  pyrometer  trace  is  shown  in  Figure  13.  Note  the  plateau  corresponding  to 
the  combined  model  surface  emission  and  radiant  source  reflected  components;  the  difference  between 
"peak-to-val ley"  levels  represents  the  reflected  component. 

3.2.4  Ablation  Mod e  1_  Design 

The  principal  thermal  consideration  of  the  model  design  was  to  reduce  heat  conduction  through 
the  model  material  and  radiation  losses  from  the  cylindrical  surface  and  thus  maintain  the  stagna¬ 
tion  area  temperature  at  the  maximum  possible  level.  This  criterion  had  to  be  met  while  maintaining 


Figure  12.  Pyrometer  alignment  -  top  view. 
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the  structural  integrity  of  the  model  during  injection  pressure  loading  and  the  initial  thermal 
shock  conditions.  Several  configurations  were  investigated  in  the  preliminary  system  checkout  tests 
by  exposing  proof  models  of  the  following  designs: 

1.  "Mushroom"  shaped  models  with  a  0.25-inch  radius  hemispherical  tip  and  a  small  (0.25- 
inch)  stem  running  the  length  of  the  model  to  the  holder.  The  stem  was  surrounded  by 
a  concentric  annular  cross  section  radiation  shield  of  the  same  material  as  the  model. 

2.  Hollow  thin-walled  models 

3.  Baseline  solid  design 

Various  stan  diameters  and  wall  thicknesses  were  investigated  in  designs  (1)  and  (2)  above  attempting 
to  1  jini t  the  conduction  losses  while  maintaining  structural  and  thermal  design  criterion.  Although 
the  solid  model  was  structurally  superior  and  easier  to  fabricate,  the  conduction  losses  were  sig¬ 
nificantly  higher  than  the  other  designs.  The  hollow  models  experienced  thermal  shock  problems 
(i.e.,  cracked)  and  proved  more  difficult  to  manufacture.  Design  one  was  fourd  to  be  optimum  from 
thermal  response  and  manufacturing  considerations;  therefore,  it  was  used  exclusively  throughout 
the  test  series. 

The  overall  model  size  was  limited  by  the  size  constraints  of  the  radiation/optical  system 
and  the  axial  movement  of  the  servo-controlled  model  positioner  system.  As  shown  in  Figure  6, 
imaging  mirror  focal  point  diameter  is  approximately  0.5  inch,  which  limits  the  model  diameter  to 
approximately  0.5  inch  for  uniform  heating  of  the  model.  The  focal  length  of  the  radiation  imaging 
mirror  and  the  axial  range  of  the  positioning  system  define  the  limits  of  the  model's  overall  length. 
The  model  length  must  be  such  that  there  is  sufficient  axial  movement  capability  to  actively  con¬ 
trol  the  stagnation  point  at  the  imaging  mirror  focal  point  throughout  the  test. 

The  test  models  were  machined  from  Union  Carbide  ATJ-S  graphite  ("P"  bil let).  They  were 
0.25  inch  hemispherical ly  tipped  cylinders  1.5  inches  long.  The  cylindrical  body  consisted  of 
an  inner  stem  and  a  surrounding  "sleeve."  Figure  14  is  a  sketch  of  the  model  assembly.  The  model 
was  mated  to  the  facility  water  cooled  sting  by  a  carbon  phenolic  model  holder  (MX-4926  or  FM5790). 

The  model  was  cross-pinned  to  the  holder  and  the  holder  was  in  turn  attached  to  the  sting  by  a  set 
screw. 

3.3  TEST  PLAN 

The  test  program  was  completed  in  two  distinct  phases.  The  first  phase,  cneckout,  proof  and 
calibration  tests,  was  designed  to  "debug"  the  operation  of  the  arc  heater  and  is  described  in  Section 
3.3.1.  The  second  phase  consisted  of  the  actual  sublimation  tests,  as  discussed  in  Section  3.3.2. 
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3.3.1  Checkout,  Proof  and  Calibration  Tests 

The  objective  of  this  initial  test  series  was  to  fully  characterize  the  operation  of  the 
new.y  configured  arc  heater  while  generating  the  data  necessary  to  select  the  model  design  and  fa¬ 
cilities  optimum  operating  conditions.  Since  this  arc  heater  had  never  been  operated  in  the  AEHS  fa¬ 
cility,  and  further,  the  expected  convective  levels  were  much  higher  than  those  produced  with  pre 
vious  heaters,  certain  operational  problems  had  to  be  overcome  in  the  initial  checkout  tests.  These 
problans  are  summarized  as  follows: 

1.  Verification  of  arc  starting  technique 

2.  Verification  of  power  supply  stability  on  arc  initiation  (i.e.,  no  current  surge) 

3.  Constantly  rising  test  cabin  pressure  due  to  nature  of  vacuum  system  and  cooling/purge 

flows 

The  arc  heater  was  initially  configured  in  the  longest  (46  segment)  arrangement  and  poten¬ 
tially.  the  most  difficult  starting  condition.  The  first  starting  attempt  resulted  in  a  large  cur¬ 
rent  surge  with  subsequent  failure  of  arc  heater  components  and  an  external  arc-over.  After  making 
adjustments  to  the  power  supply  controls  to  prevent  further  current  surges  of  this  magnitude,  the 
arc  heater  was  rebuilt  with  the  following  modifications: 

1.  All  potential  arc-over  sites  around  and  on  the  arc  heater  were  carefully  insulated 

2.  A  small  amount  of  the  primary  test  gas  (normally  N2,  but  argon  during  the  start  sequence) 
was  bypassed  and  injected  over  the  face  of  the  cathode  to  increase  arc  stability 

A  subsequent  checkout  test  demonstrated  the  proper  arc  operation  with  this  46  segment  con¬ 
figuration.  The  arc  was  started  in  this  configuration  through  the  use  of  a  small  argon  flow  as  the 
primary  gas  and  sufficient  power  supply  open  circuit  voltage  to  initiate  the  arc.  Once  the  arc  was 
stable,  the  required  primary  gas  (N2)  was  introduced  and  the  argon  flow  terminated.  Later  the  oxy¬ 
gen  was  introduced  to  provide  proper  simulation  of  air. 

The  heater  was  then  reconfigured  into  the  23  segment  constrictor  length  (Configuration  3). 
initial  plans  called  for  a  separate  checkout  series  followed  by  a  short  calibration  series  to  define 
four  test  conditions  to  be  used  in  the  subsequent  proof  tests.  The  proof  tests  were  designed  to 
expose  several  different  model  designs  so  that  the  optimum  model  design  would  be  selected  for  the 
primary  test  series.  Due  to  the  unforeseen  length  of  time  required  to  accomplish  each  test  point, 
the  checkout  test  series  was  modified  significantly.  For  each  arc  heater  configuration,  checkout, 
calibration  and  proof  test  were  to  be  combined  into  a  short  two  or  three  test  sequence.  In  this 
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manner,  problems  encountered  later  in  the  checkout/calibration/proof  test  sequence  would  not  pre¬ 
vent  the  successful  completion  of  at  least  part  of  the  sublimation  tests. 

With  the  arc  configured  with  23  segments  (Configuration  3)  the  test  cabin  pressure-rise  rates 
were  studied  while  calibration  data  were  taken  for  the  first  proof  test.  Six  runs  were  made  to 
study  ways  of  decreasing  the  test  cabin  pressure-rise  rate  and  its  effect  on  the  measured  heating 
rate  and  stagnation  pressure.  The  final  two  runs  in  this  sequence  provided  stream  profiles  used 
for  final  quantification  of  the  effect  of  the  pressure  variation  on  stream  properties. 

At  this  point  In  the  preliminary  test  series  the  radiant  source  was  calibrated  and  the  first 
proof  tests  run.  A  total  of  four  models  were  tested  (Runs  18  and  20)  two  per  run.  The  first  of 
these  became  a  convection-only  test  when  the  radiant  source  failed  to  operate  throughout  the  entire 
test.  The  second  attempt  was  successful  with  good  pyrometer  and  film  data  being  obtained. 

With  the  completion  of  the  first  proof  tests,  facility  scheduling  requirements  forced  a  halt 
to  any  further  tests.  These  few  tests  however,  did  provide  the  following  information. 

1.  The  variation  of  the  test  cabin  pressure  results  in  variation  of  test  conditions  as  the 
shock/ expansion  structure  of  the  stream  changes  (applies  to  supersonic  mode). 

2.  The  Thermogage  pyrometer  tracked  the  chopper  very  well.  Alignment  was  shown  to  be  ex¬ 
tremely  important  due  to  the  erratic  behavior  of  some  of  the  temperature/time  history 

plots. 

3.  All  model  designs  were  structurally  adequate,  although  the  hollow  models  showed  a  ten¬ 
dency  to  crack. 

4.  The  model  holder  was  redesigned  to  prevent  damage  to  the  facility  sting.  The  flare  anale 
was  increased  and  a  boron  nitride  shield  was  added  to  facilitate  removal  of  model  from 
the  sting. 

3.3.2  Sublimation  Tests 

With  only  one  of  the  six  arc  heater  configurations  fully  calibrated,  as  described  in  the  pre¬ 
vious  section,  scheduling  restrictions  imposed  by  the  AEHS  facility  resulted  in  modification  to  the 
sublimation  test  sequence.  It  became  necessary  to  conduct  simultaneous  calibration/proof  tests 
and  primary  sublimation  tests  at  a  specific  condition.  This  section  includes  (1)  a  description  of 
the  testing  sequence  employed  and  (2)  review  of  specific  problems  and  their  solutions  encountered 
during  this  second  entry  into  the  AEHS  facility.  Also  presented  is  a  summary  and  data  tabulation 
for  each  of  the  ablation  test  conditions  run. 
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The  sequential  calibration/proof  test  procedure  resulted  in  several  of  the  early  runs  (up 
through  run  47)  being  far  from  optimum  with  respect  to  the  test  objects  of: 

Maximum  accuracy  of  surface  temperature  data  which  implies  the  highest  possible  surface 
temperature  measurements. 

this  problem  did  dot  affect  the  later  rent  becaete  of  procedural  changed  in  pyrometer  alignment 


techniques. 

After  reconfiguring  the  arc  for  a  specific  test  condition  the  procedure  followed  involved 
(,)  verifying  the  arc  starting  techniques  and  (2)  reviewing  calibration  runs  to  determine  the 

test  condition  parameters. 

The  typical  test  sequence,  after  preliminary  calibration  data  had  characterized  the  test 
condition,  was  as  follows: 

A.  Radiant  source  started  and  set  to  the  desired  source  current,  -ithout  the  chopper  running, 
radiant  heating  rate  data  »ere  taken  with  the  steady-state  calorimeter  and  the  blackened 

slug  calorimeter- 


li.  Douser  was  closed  and  the  chopper  started. 

C.  Argon  start  flow  initiated  and  the  arc  heater  started.  With  arc  running  in  a  stable 
mode,  N2  and  02  were  injected  and  the  argon  flow  terminated. 

0.  After  the  predetermined  arc  current  had  been  achieved  and  a  stable  arc  established,  the 
convective  conditions  were  measured  with  the  steady-state  calorimeter  and  pitot  probe. 

E.  The  douser  was  opened  and  the  high  speed  camera  started  prior  to  the  model  arriving  on 
the  stream  centerline.  The  ablation  model  was  positioned.  Ablation  test  duration  was 
automatically  controlled  to  a  preset  value,  with  douser  closure  and  removal  of  the  model 

from  the  test  stream. 

F.  While  the  arc  heater  continued  to  run  at  the  test  condition,  post-test  data  were  taken 
with  the  steady-state  calorimeter  and  the  pitot  probe. 

G.  The  arc  heater  was  turned  off,  the  chopper  stopped,  the  douser  opened  and  a  final  radiant 
measurement  taken  with  the  steady-state  calorimeter. 


H.  Radiant  source  turned  off. 

I.  Test  chamber  vented  to  atmospheric  pressure.  Ablation  model  repositioned  and  post-test 
alignment  check  of  pyrometer  made  and  noted.  Model  removed  from  sting. 
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Model  *13 


Pre-  and  post-test  data  were  necessary  to  completely  define  the  test  environment  in  the  chang¬ 
ing  test  cabin  pressure  situation.  While  considerable  work  was  done  duirng  the  checkout/calibration 
phase  to  minimize  the  pressure  rise  during  the  nominal  2-minute  test  sequence,  there  remained  a 
significant  effect  on  the  convective  environment.  This  was  due  mainly  to  the  change  in  the  test 
stream  shock/expansion  Wave  structure  in  the  istper sonic  case,  i.e.,  the  test  rhombus  continually 
changed.  Test  cabin  pressure  variations  in  the  subsonic  case  directly  affected  the  model's  stag¬ 
nation  pressure. 

Operation  of  the  radiation  source  was  steady  as  monitored  by  the  controls  for  the  argon  arc 
source.  Pre-  and  post-test  radiation  calibration  data,  however  showed  a  consistent  decay  in  the 
incident  radiation  due  to  (1)  decay  in  the  reflectance  of  the  imaging  mirror  due  to  soot  deposition 
and  (2)  a  general  degradation  in  the  transparency  of  the  static  gas  environment  due  primarily  to 
soot  arid  N02  formation  in  the  test  chamber.  The  incident  radiation  data  reported  reflect  this  deg¬ 
radation. 

The  0.1  atm,  supersonic  sublimation  tests  are  summarized  in  Table  4.  The  heating  rate  and 
stagnation  pressure  data  are  the  average  of  the  pre-  and  post-test  measurements.  The  centerline 
enthalpy,  Hq^  ,  is  also  an  average  based  on  the  pre-  and  post-test  calorimeter  measurements,  as  de¬ 
termined  by  a  procedure  described  in  Section  3. 4. 2.1.  These  comments  apply  to  the  summary  tables 
presented  in  this  section.  The  complete  set  of  data  are  shown  in  Appendix  B.  Test  22  was  a  final 
arc  heater  calibration  of  this  condition,  with  radiant  source  calibration  being  accomplished  in  tests 
23  and  24.  Initial  sublimation  data  were  obtained  from  tests  25,  26  and  28.  Subsequent  analysis  of 
these  data  revealed  that  the  surface  temperature  measured  was  significantly  lower  than  expected.  Two 
possible  causes  were  identified: 

1.  Misalignment  of  pyrometer 

2.  High  conduction  losses  within  the  model  (insufficient  incident  radiative  heating  rates) 

The  radiant  source  current  was  increased  to  the  1300  to  1500  amp  level  to  maximize  the  radiant 
energy  flux  to  the  model  surface.  In  addition  great  care  was  taken  in  aligning  the  pyrometer.  Sub¬ 
sequently  tests  58,  61  ,  62,  63  and  64  were  run.  These  tests  provided  excellent  ablation  data  points. 

The  0.1  atm,  high  velocity,  subsonic  sublimation  tests  are  summarized  in  Table  5.  Since  no 
checkout/calibration  tests  had  been  run  at  this  test  condition,  test  32  (no  data  shown)  was  the  ini¬ 
tial  checkout  test  for  the  arc  heater.  Tests  33  to  37  were  calibration  runs,  with  various  arc  heater 
flowrates  and  arc  currents  being  used  to  define  the  optimum  arc  combination.  Since  the  test  stream 
was  subsonic,  changes  in  the  test  cabin  pressure  were  directly  translated  as  changes  in  the  model 
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stagnation  pressure,  which  accounts  for  the  large  test  to  test  pressure  variation.  Initial  sub¬ 
limation  data  were  obtained  from  tests  38  and  39.  For  the  reasons  outlined  above,  these  tests  were 
rerun  as  tests  48,  49  and  50. 

The  0.1  atm,  low  velocity,  subsonic  sublimation  tests  are  summarized  in  Table  6.  The  cali¬ 
bration  and  initial  sublimation  tests  were  tests  43  to  46.  Final  sublimation  data  were  obtained 
from  tests  51  to  54.  Ablation  data  from  these  tests  are  presented  in  Section  3.4.3. 

The  1.0  atm,  supersonic  sublimation  tests  are  summarized  in  Table  7.  At  this  condition,  the 
test  cabin  pressure  rise  had  the  least  effect  on  the  convective  test  stream  due  to  the  relatively 
high  pressure  (70  to  80  torr)  required  for  matched  nozzle  flow.  This  allowed  for  a  lower  initial 
cabin  pressure,  maintaining  operation  in  the  underexpanded  mode  throughout  the  test.  The  result  was 
a  test  rhombus  which  varied  only  slightly  over  the  test  period.  Test  69  provided  arc  heater  cali¬ 
bration  data  while  final  sublimation  data  was  obtained  from  tests  70,  71  ,  72  and  73.  A  total  of  4 
usable  data  points  were  obtained  at  this  condition. 

The  final  sublimation  tests,  1.0  atm,  low  velocity,  subsonic,  are  summarized  in  Table  8.  For 
this  condition,  th  >  test  cabin  was  vented  to  the  atmosphere,  thus  maintaining  the  stagnation  pres¬ 
sure  at  1.0  atm.  The  high  cabin  pressure  produced  an  unexpected  side  effect.  The  production  of 
various  nitric  oxide  compounds  ( N0x )  was  very  high  and  appeared  to  be  predominately  NO.,  from  its 
reddish-brown  color.  It  is  felt  this  was  a  result  of  the  high  oxygen  content  of  the  cabin  at  this 
condition  reacting  with  the  high  temperature  nitrogen/oxygen  test  stream.  Since  this  operating  mode 
was  unique,  there  was  no  means  of  evacuating  the  cabin  during  the  test  at  a  rate  sufficient  to  pre¬ 
vent  N0X  formation.  Visual  observation  of  tests  76,  77  and  78,  which  were  run  with  arc  heated  oxy¬ 
gen  and  nitrogen,  revealed  a  significant  reduction  in  the  radiant  levels  due  to  the  N0x  buildup. 

Pre-  and  post-test  measurements  of  the  radiant  levels  (see  Tahle  8)  showed  a  40  to  50  percent 
reduction.  Tests  79  and  80  were  an  atteupt  to  reduce  the  concentration  of  N0X  by  removing  the  arc 
heated  oxygen.  These  two  tests  were  made  using  arc  heated  nitrogen.  This  modification  resulted  in 
a  radiant  reduction  of  33  to  35  percent,  indicating  that  arc  heated  nitrogen  was  still  enabling  the 

N0X  reactions  to  proceed  at  a  significant  rate.  It  should  be  noted  that  the  N0X  buildup  in  the 

cabin  not  only  reduced  the  radiant  levels  incident  on  the  model  but  also  perturbed  the  pyrometer 
data  (see  Section  3.4.1).  For  this  reason,  the  final  test  was  run  with  cold  argon  flowing  over  the 
mode!  at  a  low  subsonic  velocity  of  about  2  ft/sec;  the  only  energy  input  to  the  ablation  model  was 

radiant.  The  intention  of  the  argon  flow  was  to  remove  ablation  products  (i.e.,  soot  formed  by  con¬ 

densing  carbon  vapor  species)  from  the  vicinity  of  the  ablating  surface  which  block  some  of  the 
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emitted  radiation.  In  this  test,  (#81)  the  argon  flow  elocity  was  apparently  too  low  to  be  effec¬ 
tive  and  a  plume  or  cloud  of  carbon  soot  was  observed  at  the  model  surface.  The  effect  of  this 
plume  is  discussed  in  Section  3.4. 

In  summary,  sublimation  data  points  were  obtained  at  two  pressure  levels  and  three  distinct 
flow  velocities  for  a  total  of  five  separate  convective  conditions. 

3.4  TEST  RESULTS 

The  results  of  the  graphite  vaporization  study  conducted  in  the  NASA  Ames  AEHS  facility  are 
presented  in  this  section  along  with  a  discussion  of  the  techniques  used  to  reduce  the  data.  The 
test  matrix  with  a  tabulation  of  measured  test  conditions  is  presented  in  Section  3.4.1.  The  pro¬ 
cedures  used  to  characterize  the  convective  environment,  verify  the  pyrometer  data,  and  determine 
the  surface  recession  rates  are  discussed  in  Section  3.4.2.  Uncertainties  associated  with  each 
type  of  data  are  analyzed.  The  ablation  response  of  each  test  model  is  presented  in  Section  3.4.3. 

3.4.1  Test  Matrix  and  Test  Condi  cions 

Table  9  is  the  test  matrix  and  tabulation  of  the  average  stream  and  model  surface  test  con¬ 
ditions  for  all  ATJ-S  ablation  model  exposures  in  this  test  series.  The  initial  proof  model  ex¬ 
posures  and  calibration  runs  are  not  presented,  but  can  be  found  in  Appendix  B  along  with  actual 
stream  measurements  before  and  after  model  exposures.  The  complete  tabulation  in  Appendix  B  also 
describes  the  assumptions  required  for  those  data  that  were  unavailable  due  to  data  acquisition  or 
probe  malfunction. 

3.4.2  Data  Reduction  Techniques 

The  data  reduction  procedures,  assumptions  and  uncertainties  involved  in  convective  stream 
characterization,  model  surface  temperature  measurements,  and  recession  rates  are  discussed  in  the 
following  subsections. 

3. 4. 2.1  Convective  Environment  Characterization 

As  mentioned  earlier  in  Section  3.2.3,  the  arc  heater  cooling  water  temperature  rise  and  mass 
flow  rate  were  not  measured  in  this  facility  and,  therefore,  the  arc  heater  efficiency  and  bulk 
enthalpy  were  not  known.  In  addition,  the  low  heating  rates  encountered  in  the  subsonic  portions 
of  this  test  series,  coupled  with  the  necessity  to  use  a  0-5000  Btu/ft2  sec  steady-state  calorimeter 
(see  Section  3.2.3)  resulted  in  cold-wall  convective  heating  measurements  corresponding  to  less 
than  1  percent  of  the  full  range  output  of  the  instrument.  These  two  factors  affected  the  charac¬ 
terization  of  the  convective  environment  at  the  subsonic  conditions.  Under  both  the  supersonic  and 
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52 

29 
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53 

20 

29.35 

2890 

50.5 
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33 
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! The  copper  steady-state  calorimeter  was  not  blackened  -  the  q  data  were  adjusted  by  an  emissivity 
of  0.83  based  on  blackened  transient  calorimeter  data. 

*Based  on  Fay-Riddell  theory  of  stagnation  point  convective  heat  transfer  to  a  5/16"  calorimeter. 


subsonic  flow  conditions  the  arc  heater  chamber  pressure  (Pc)  was  accurately  measured  by  a  pressure 
transducer.  The  stream  stagnation  pressure  (Pt)  was  measured  with  a  pitot  probe  inserted  into  the 
stream  prior  to  and  immediately  following  the  test  model  exposure.  The  determination  of  the  con¬ 
vective  flow  centerline  enthalpy  was  based  on  the  measured  cold  wall  heat  flux  but  was  calculated 
differently  for  the  subsonic  and  supersonic  cases.  These  two  procedures  are  described  below. 


Supersonic  Condition 


ror  tne  two  supersonic  test  conditions  (arc  heater  coTif  iyo,  atioTiS  3  LTid  €),  the  tbfiterllfit 
enthalpy  was  calculated  using  the  Fay-Riddel  1  theory  for  cold  wall  stagnation  point  convective  heat 
transfer  to  a  sphere  and  the  measured  cold  wall  heat  flux.  The  heating  rate  (q)  measured  by  the 
s..eaoy-stcii.fc  calur  imetfci  was  iii  toe  mid  to  uppeT  ranyt  of  tiit  tiiStttnirent .  Thesv  Jata  along  with 
the  measured  stagnation  pressure  (P^)  were  use(1  t0  compute  the  centerline  enthalpy  (Hq  )  from  the 
following  relation  (Reference  43): 


q  »  0.0462  •''P^/Rn  ^o 

where:  R^  =  nose  radius  of  the  steady-state  calorimeter  (ft) 
Pt2  =  stagnation  pressure  (atm) 
q  =  heating  rate  (Btu/ft2-sec ) 


Hq  =  stream  enthalpy  (Btu/lbm) 

The  measured  q  data  at  these  test  conditions  are  considered  accurate  and  were  used  to  determine  the 
effective  centerline  enthalpy  for  each  supersonic  run. 


Subsonic  Conditions 

The  subsonic  stream  enthalpy  was  calculated  using  an  iteration  technique  which  required  an 
arc  heater  efficiency  to  be  assumed  and  the  bulk  enthalpy  (HQ)  calculated  usling  the  known  input 
power  and  gas  mass  flow  rate.  This  initial  enthalpy  estimate  along  with  the  measured  chamber 
pressure  (Pc)  was  used  to  define  the  arc  chamber  gas  state.  At  these  subsonic  flow  conditions, 
due  to  uncertainties  in  the  measured  data,  the  data  from  all  the  runs  at  a  particular  subsonic  con¬ 
dition  were  used  to  define  an  average  Pc,  HQ  and  cold  wall  heat  rate  (q ) .  The  subsonic  heat  transfer 
coefficient  to  an  axisymmetric  blunt  body  (Reference  44)  was  computed  and  compared  to  the  experi¬ 
mentally  determined  value  based  on  the  average  measured  q  and  initial  effective  H  based  on  an 
assumed  efficiency.  Various  efficiencies  were  assumed  in  this  iteration  until  closure  was  achieved 


and  the  theoretical  and  experimental  heat  transfer  coefficients  agreed.  The  following  relations 
were  used  in  the  theoretical  evaluation  /if  the  heat  transfer  coefficient: 


I  I 

Nux  -0.76  Rex  lj  I>r°  -=  h^x 


These  relations  assume  (1)  incompressible  flow,  (2)  no  centerline  enthalpy  peaking,  and  (3)  a  fully 
filled  nozzle. 

As  previously  mentioned,  the  low  heating  rates  of  the  subsonic  test  conditions  and  the  rela¬ 
tively  high  range  of  the  steady-state  calorimeter  caused  substantial  uncertainty  in  the  subsonic 
heating  data.  Therefore,  when  the  average  q  for  runs  of  similar  conditions  at  each  subsonic  test 
point  were  computed,  obviously  errorous  data  were  ignored.  The  similarity  in  material  response, 
power  input,  stream  conditions  and  available  convective  q  data  from  the  transient  slug  calorimeter 
were  used  to  define  an  average  q  for  each  subsonic  condition.  Arc  heater  Configuration  1  is  an 
example  of  the  selection  processes  used  to  average  the  subsonic  q  data.  Of  the  four  runs  at  this 
condition  (runs  51  through  54),  only  three  produced  a  recorded  q  data.  Only  two  of  these  were  con¬ 
sidered  valid  (runs  52  and  53).  The  low  q  data  from  run  54  were  ignored  (see  Table  9).  The  low 
millivolt  output  resulted  in  a  barely  discernible  deflection  on  the  CEC  recording  system.  However, 
the  surface  recession  rates  of  the  four  runs  were  similar.  Since  the  repeatability  of  the  measured 
q  (50.5  Btu/ft-sec)  was  evident  throughout  the  four  calorimeter  exposure  of  runs  52  and  53,  this 
value  was  assumed  as  the  average  q  for  Configuration  1.  Since  the  sensitivity  of  the  calorimeter 
was  low  at  these  test  conditions,  the  estimated  accuracy  of  the  subsonic  q  data  due  to  deflection 
measurement  uncertainties  alone  is  *20  percent.  Although  these  uncertainties  are  large  their  effect 
on  the  data  interpretation  isn't  too  significant  as  shown  in  Section  4. 
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The  last  series  of  runs  at  the  1.0  atm  subsonic  configuration  produced  equally  suspect 
steady- state  calorimeter  data.  For  these  runs  the  average  convective  q  measured  by  the  transient 
slug  calorimeter  was  also  used  to  estimate  the  stream  enthalpy  and  subsonic  heat  and  mass  transfer 

rate. 

3.4. 2.2  Pyrometer  Data  Reduction 

Several  potential  problem  areas  that  affect  the  optical  pyrometer  surface  temperature  mea¬ 
surements  were  investigated  to  determine  their  impact  on  the  reported  data.  Potential  problems 
include  alignment,  gas  cap  radiation,  ablation  products  in  the  boundary  layer,  and  the  occlusion  of 
tho  reflected  component.  Each  of  these  areas  is  discussed  in  this  section. 

A1  i  gnment 

During  the  initial  runs  of  the  test  program  pyrometer  alignment  became  a  series  problem  and 
resulted  in  the  low  temperature  data  reported  in  runs  38  through  45.  Initially  the  target  was  po¬ 
sitioned  off  of  the  stagnation  point  along  the  45°  ray  in  the  horizontal  plane  (see  Figure  12). 
Thermal  and  pressurization  cycling  of  the  models,  as  well  as  slignt  variations  in  model  length  and 
positioning  differences,  resulted  in  rather  gross  uncertainties  in  the  pyrometer  target  location  on 
the  model  surface.  When  this  fluctuation  in  temperature  became  apparent,  greater  care  was  taken  in 
verifying  the  pyrometer  alignment.  The  target  location  was  checked  and  noted  prior  to  each  run. 
After  the  model  exposure  the  model  was  reinserted  and  the  target  location  checked  again.  The  loca¬ 
tion  with  respect  to  the  maximum  recession  area  was  noted.  These  steps  resulted  in  highly  consis¬ 
tent  temperature  levels  for  similar  test  conditions  for  runs  48  through  the  completion  of  the  test 
series. 


Gas  Cap  and  Ablation  Products  Radiation 

A  potential  source  of  error  in  pyrometer  data  of  ablating  graphite  surfaces  is  the  emission 
and/or  absorption  of  the  high  temperature  gas  cap  between  the  normal  shock  and  the  model.  Another 
concern  is  the  radiative  properties  of  the  ablation  products  within  the  boundary  layer.  The  two 
effects  are  shown  to  be  negligible  by  Rindal,  et  al .  (Reference  45)  for  carbon  sublimation  in  air 
under  environmental  conditions  similar  to  these. 

in  Reference  45  gas  cap  radiation  was  calculated  for  stagnation  pressures  (Pt2)  of  0.05  atm, 

0.5  atm,  and  100  atm  assuming  air  enthalpies  of  30,000  Btu/lbm  at  the  lower  pressures  and  3000  Btu/lbm 
at  the  high  pressure  condition. 
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In  calculating  the  radiation  from  the  gas  cap,  the  gas  cap  was  assumed  to  be  optically  thin 
with  the  radiation  being  a  combination  of  both  continuum  and  band  (line)  radiation  The  thickness 
of  the  gas  cap,  the  shock  stand  off  distance,  was  calculated  from: 

P 

4  -  0.78  —  Ru 
po  N 

where  R^  =  nose  radius,  0.5  inch  for  the  AFML  study  note  that  the  nose  radius  for  all  models  of 
the  Aerotherm  study  reported  herein  was  0.25  inch 

P„  =  free  stream  density 

PQ  =  stagnation  point  density 

The  radiation  from  the  ablation  products  was  studied  for  only  one  condition  P.  =  0.5  atm  and 
Hq  -  30,000  Btu/lbm.  The  effects  of  both  the  ablation  products  and  the  gas  cap  are  summarized  in 
Figure  15.  The  emissive  power  versus  wavelength  from  the  gas  cap  for  two  conditions  and  that  from 
the  ablation  products  at  one  condition  are  compared  to  the  emissive  power  of  a  bla. kbody  surface  at 
various  temperatures.  The  hand  radiation  component  is  superimposed  upon  the  continuum  radiation; 
the  total  radiation  power  being  equal  to  the  sum  of  these  two  levels. 

As  shown  in  Figure  15,  for  wavelengths  greater  than  0.2p  the  bard  contribution  wa-;  negligible 
compared  to  the  continuum  contribution  for  both  the  high  and  low  pressure  conditions.  Furthermore, 
at  the  100  atm  condition,  above  0.5g,  the  gas  cap  effects  were  less  than  or  equal  to  ’c. 6  percent  of 
the  emissive  power  of  a  blackbodv  at  7000“R  (the  nominal  range  of  temperatures  reported  in  this  in¬ 
vestigation).  The  test  pressures  of  this  study  were  in  the  range  (0.1  to  1.0  atm)  which  is  close  to 
the  0.5  atm  curve  of  Figure  15,  where  it  can  be  seen  that  the  gas  cap  effects  are  even  less  pronounced. 

The  ablation  products  emissive  power  contribution  is  less  than  1  percent  above  0.6p  relative 
to  a  7000°R  blackbody .  Radiation  from  the  freestream  gas  upstream  of  the  shock  was  also  found  to  be 
negligible  in  this  study. 

Reflected  Component 

As  discussed  previously  in  Section  3.2.1,  the  combined  radiation  facility  utilized  a  large 
diameter  rotating  chopper  to  occlude  the  incident  radiation  from  the  radiation  source  and  thus  elimi¬ 
nate  the  model  surface  reflectance  component  from  the  pyrometer  measurement.  This  technique  is  only 
successful  if  the  pyrometer  and  associated  recording  apparatus  response  time  is  sufficiently  fast 
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igure  15.  Radiation  contributions  of  surface  gas  cap  and  ablation 
products  from  AFML-TR-70-133. 


to  monitor  the  pyrometer  output  during  the  short  blockage  time.  The  chopper  blades  occlude  the  inci¬ 
dent  radiation  for  3  msec  with  an  interval  of  about  50  to  60  msec  between  each  blockage.  The  10 
psec  response  of  the  pyrometer  is  more  than  required  for  this  operation.  However  the  response  time 
of  the  recording  equipment  is  such  that  the  effective  response  time  of  the  system  is  close  to  3  msec. 

The  response  time  of  the  pyrometer  and  recording  system  was  demonstrated  by  exposing  the 
water  cooled  steady-state  calorimeter  to  the  incident  radiation  environment  while  the  chopper  was 
operating.  The  pyrometer  output  dropped  to  the  cold  wall  value  of  the  calorimeter  during  the  time 
when  the  chopper  blades  blocked  the  incident  radiation.  This  response  is  shown  In  Figure  16. 

34.2.3  Recession  Data 

Model  centerline  recession  rates  were  generated  by  reducing  the  high  speed  16mm  film  data  of 
each  model  exposure.  An  L&W  photo-optical  data  analyzer  was  used  to  project  the  model  images.  The 
location  of  the  model  centerline  with  respect  to  a  fixed  model  reference  line  was  noted  at  several 
times  throughout  the  test.  The  times  for  each  location  were  computed  by  defining  an  average  framing 
rate  from  the  projector  frame  counter  and  the  overall  test  time  from  the  CEC  oscillograph  trace.  The 
recession  to  each  profile  was  computed  by  measuring  the  projected  images  with  respect  to  a  fixed 
reference  and  converting  to  actual  dimensions  through  the  use  of  a  scale  factor.  The  reference  used 
was  the  interface  between  the  primary  sphere/cylinder  portion  of  the  model  and  the  aft  cylindrical 
radiation  shield.  The  scale  factor  was  determined  from  an  average  of  the  projected  model  diameter 
over  several  frames  near  time  0  (after  the  model  reached  test  stream  centerline)  and  the  known 
pretest  diameter.  The  scaling  factor  was  determined  early  in  time  to  reduce  the  effects  of  sidewall 
ablation  and  thermal  expansion.  With  these  phenomena  in  mind,  noting  that  precise  measurements 
of  some  projected  images  were  difficult  due  to  film  clarity  and  contrast,  the  film  data  length  mea¬ 
surements  are  considered  accurate  to  within  about  *10  percent. 

The  data  from  the  film  reduction  procedure  were  plotted  as  recession  versus  time.  A  straight 
line  was  hand-fitted  through  the  data  over  the  steady-state  portion  of  the  run.  The  slope  was  com¬ 
puted  and  the  results  are  presented  in  the  following  subsection. 

For  several  runs  the  exposure  settings  and  filtering  resulted  in  poor  contrast  such  that  the 
model  reference  line  could  not  be  clearly  identified  during  the  intermediate  portion  of  the  test. 

For  these  cases  an  approximate  recession  rate  was  evaluated  by  the  ratio  of  the  overall  length  change 
to  the  total  test  time  (As/A0).  The  length  change  was  the  difference  between  the  initial  and  final 
lengths  as  determined  from  the  film  data.  Initial  and  final  model  dimensions  were  available  from 
the  poor  quality  film  data  since  the  reference  line  was  generally  visible  as  the  model  was  moving 
into  and  out  of  the  test  stream.  At  these  times  the  model  was  outside  of  the  high  intensity  focal 
point  of  the  radiation  source  imaging  mirror  such  that  the  photographic  contrast  was  better. 
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The  As/A0  values  as  determined  from  the  film  data  were  considered  more  accurate  than  As/A6 
values  determined  from  pre-  and  post-test  length  measurements  of  the  test  samples  for  the  following 
reasons : 

1.  The  pre-  and  post-test  measurements  were  made  on  the  overall  model  length,  including  the 
model  holder;  thermal  cycling  may  have  effected  the  model/model  holder  junction 

2.  Some  ablation  occured  after  removal  from  the  test  stream 

3.  Handling,  after  the  exposure  and  before  post-test  measurements  were  taken,  may  have 
damaged  the  surface  and  affected  the  As  measurement 

For  those  runs  where  no  film  data  are  available  the  s  data  were  estimated  from  other  runs  at 
similar  conditions.  A  correlation  between  s  as  measured  by  film  reduction  techniques  and  As/AO  as 
determined  by  pre-  and  post-test  length  measurements  was  used  to  estimate  s  from  the  As/A0  value  of 
the  effected  model. 

For  the  subsonic  runs  the  model  ablation  was  slightly  asymmetric.  This  is  evident  from  the 
photographs  in  Figure  17  of  four  ablation  models  tested  at  the  four  primary  test  conditions.  The 
area  off  of  the  model  centerline  along  the  45°  ray  receeded  more  than  the  stagnation  area.  This 
high  recession  area  was  the  viewing  location  of  the  optical  pyrometer  but  was  not  viewed  by  the  high 
speed  motion  picture  camera  (see  Figure  12,  Section  3.2.3).  Therefore,  the  centerline  recession 
rates  reduced  from  the  film  data  were  not  indicative  of  the  higher  recession  rates  in  the  location 
of  the  surface  temperature  measurements.  To  investigate  this  disparity,  the  recession  rate  at  the 
45°  location  was  estimated  from  the  ratio  of  recession  at  the  45°  ray  ( As^g )  to  centerline  recession 


The  exposed  test  models  which  exhibited  significant  asymmetries  were  projected  using  an  optical  com- 
paritor  (10X  power)  and  the  image  traced  in  a  plane  corresponding  to  the  horizontal  plane  when  in¬ 
stalled  in  the  test  facility.  This  plane  showed  the  asymetric  ablation  shape  at  the  45°  location. 
In  this  plane  the  AsQ  occuring  at  the  stagnation  point  was  approximately  equal  to  the  Asq^  as  viewed 
by  the  camera.  An  unexposed  model  was  then  projected  and  the  images  superimposed.  The  overall  re¬ 
cession  at  the  centerline  and  at  the  45°  position  were  measured  and  converted  to  actual  dimensions 
using  the  appropriate  scaling  factor.  Figure  18  shows  an  example  of  this  estimation  procedure. 

For  the  subsonic  0.1  atm  conditions,  runs  38  through  54,  the  average  uncertainty  in  the 
reported  s  data  and  hence  the  blowing  parameter  8^  was  approximately  90  percent.  It  is  important 
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Post-test  ablation  models  from  the  four  primary  test  conditions 


to  note,  however  that  despite  the  rather  large  uncertainty  in  s,  the  effect  on  the  corresponding  B' 
value  is  not  significant  since  B'  values  are  all  B'  -•  10  at  these  subsonic  test  conditions. 

3.4.3  A b  1  a tion  R esults 

The  results  of  the  measured  and  "adjusted"  data  described  in  the  previous  sections  are  re¬ 
ported  in  this  section  in  terms  of  the  ablation  response  of  the  material.  Table  10  is  a  listing  of 
all  successful  model  exposure  runs  and  the  measured  or  estimated  centerline  recession  rates  (s), 

the  heat  transfer  coefficient  (p,uCH  ),  and  the  mass  flux  (m).  The  maximum  surface  temperature, 

e  e  no 

without  the  reflected  components,  and  the  calculated  blowing  parameter  (B^  =  m/p^C^)  are  also 
tabulated.  These  data  were  derived  from  the  test  conditions  reported  in  Table  9  of  Section  3.4.1. 


SECTION  4 


SUBLIMATION  DATA  ANALYSIS 

This  section  reviews  analyses  of  the  ATJ-S  ablation  test  results  presented  in  Section  3.  In 
this  section  the  ablation  data  from  the  AEHS  test  series  are  compared  with  various  thermochemical 
models  to  illustrate: 

1.  Their  agreement  with  the  JANAF  carbon  vapor  thermochemical  property  data,  and 

\  Their  correlation  with  both  equilibrium,  diffusion-controlled  and  kinetically  controlled 
models. 

The  ablation  data  generated  on  this  program  are  unique  since  tests  were  run  at  two  pressures 
(  0.1  and  1.0  atmospheres)  and  two  different  flow  conditions  (both  subsonic  and  supersonic)  at 
each  pressure.  As  discussed  in  Section  3,  the  variation  of  the  convective  flow  achieved  with  the 
modified  constrictor  arc  heater  resulted  in  about  a  factor  of  AO  change  in  the  heat  and  mass  transfer 
rate.  This  large  variation  in  pgue  Cm  is  sufficient  to  cause  about  a  200°R  -  300°R  shift  to  higher 
surface  temperatures  if  sublimation  is  rate  controlled.  Section  4.1  compares  the  ablation  data 
from  the  low  subsonic  (i.e.,  equilibrium-diffusion  controlled)  test  conditions  with  various  vapor 
pross lire  laws.  Then,  in  Section  4.2,  the  ablation  data  are  compared  with  the  equilibrium  and  rate 
controlled  thermochemical  ablation  predictions  based  on  JANAF  thermochemistry. 

4.1  DATA  COMPARISON  WITH  EQUILIBRIUM  CARBON  VAPOR  PHASE  MODELS 

The  phase  diagram  of  carbon  continues  to  be  uncertain  as  discussed  in  Section  2.  The  un¬ 
certainty  lies  principally  in: 

1.  The  thermochemical  properties  of  the  carbon  vapor  species  above  C2  (i.e.,  C3~C5)  and 

2.  The  triple  point  state. 

Results  of  the  subsonic  flow  ablation  tests  are  appropriate  for  establishing  the  correct  carbon 
vapor  species  thermochemical  model  because  of  their  low  heat  and  mass  transfer  rates  (which  re¬ 
quire  that  the  ablation  response  be  equilibrium  diffusion  controlled).  In  addition,  these  results 
can  be  compared  with  previously  reported  carbon  melt  data. 


One  means  of  establishing  the  correct  carbon  vapor  thermochemical  properties  would  be  to 
experimentally  construct  the  B'/T  curve  at  a  fixed  pressure  under  a  low  convective  flow  condition. 

The  B1  curvp  would  be  experimentally  developed  by  varying  the  incident  radiative  heating  so  as  to 
generate  carbon  ablation  data  in  the  temperature  regime  of  interest  as  previously  shown  in  Figure  3. 

Another  means  of  establishing  the  correct  carbon  vapor  thermochemistry  would  be  to  estab¬ 
lish  the  asymptotic  temperature  at  a  given  pressure  (i.e.,  high  B'  values).  At  sufficiently  high 
(S'  values  the  mass  fraction  of  carbon  vapor  at  the  surface  approaches  unity,  as  shown  by  the  ex¬ 
pression 

K^w  =  developed  in  Section  2.2. 

Thus,  for  B'  >  10,  the  carbon  vapor  partial  pressure  is  more  than  90  percent  of  the  static  pressure 
and  the  surface  state  is  controlled  by  molecular  diffusion  and  chemical  equilibrium. 

The  second  method  was  used  in  this  program  to  establish  the  eaui librium  carbon  vapor  state 
due  to  the  relatively  limited  number  of  models  and  tests  available.  At  both  pressures  (  0.1  and 

1.0  atmospheres) ,  I  he  peak  incident  radiative  flux  was  used  to  maximize  mass  flux  from  the  model 
(i.e.,  B')-  Reference  to  Table  10  shows  that  at  0.1  atmospheres  the  peak  B'  values  are  greater 
than  1 0 3 .  It  should  be  noted  that  the  blowing  correction  used  to  arrive  at  these  high  B'  values  is 
uncertain  at  these  high  blowing  conditions.  This  uncertainty  is  not  significant  however,  for  our 
purposes  since  the  accuracy  of  the  B'  value  is  not  critical  as  long  as  the  B'  value  is  greater  than 
10.  Since  the  unblown  B^  values  at  the  -  0.1  atmosphere  condition  are  above  6  it  is  obvious  the 
blown  B'  values  are  sufficiently  large  to  assume  the  carbon  vapor  is  in  equilibrium  with  the  sur¬ 
face.  Interestingly  the  B'  values  at  the  1.0  atmosphere  test  condition  are  much  lower.  This 
occurs  for  two  reasons.  first,  at  1.0  atmosphere  the  vapor  phase  equilibrium  temperature  is  600" R 
above  the  0.1  atmosphere  temperature.  The  incident  radiative  flux  was  tne  same  at  both  pressures 
(  2200  Btu/ft?  sec)  thus  the  surface  energy  balance  dictates  that  the  peak  B'  value  attainable  at 

1.0  atmosphere  must  drop  appreciably.  Secondly,  the  ablation  products  at  1  atmosphere  formed  an 
opaque  smoke  pattern  much  like  that  of  a  burning  cigarette  which  blocked  some  of  the  incident 
radiation.  In  addition  NO,,  was  formed  and  blocked  incident  radiation.  As  a  result,  the  peak  II' 
values  achieved  at  the  1  atmosphere  test  condition  are  in  the  ranqe  1  to  3.  Davy  and  Bar-Nun 
(Reference  38)  reported  a  similar  soot  or  smoke  plume  during  their  radiation  only  tests  in  the  AEHS 
facility.  In  an  attempt  to  reduce  or  remove  this  plume  on  this  series  a  low  subsonic  flow  of  arc 
heated  air  was  included  to  reduce  the  incident  and  reflected  radiation  blockage.  It  is  evident 
from  the  1  atmosphere  data  that  the  air  flow  was  insufficient  to  completely  remove  the  plume  effects 


since  the  measured  surface  temperatures  are  300°F-500°R  below  the  values  measured  under  supersonic 
flow  conditions.  In  a  further  attempt  to  alleviate  this  problem  the  last  run  (Test  81)  used  pure 
>'rt|on  as  the  convert  ive  gas  to  reduce  the  model's  m.iss  loss  (i.e.,  the  source  of  the  plume) 
and  the1  potential  for  N02  formation.  I  he  surface  temperature  measured  lor  this  run  is  about  2/b"  It 
above  the  temperatures  measured  for  those  runs  using  arc  heated  air  or  nitrogen.  In  addition  the 
B1  value  is  about  a  factor  of  4  above  those  with  the  arc  heated  air.  Therefore  it  appears  that  if 
this  technique  could  have  been  perfected  high  quality  temperature  data  would  have  been  obtained  at 
1  atmosphere  "equilibrium"  test  condition.  As  it  stands,  however,  these  surface  temperature  data 
are  felt  to  be  low. 

These  equilibrium  sublimation  data  are  compared  with  four  carbon  phase  equilibrium  predic¬ 
tions  in  F  gure  19.  The  equilibrium  data  in  the  region  of  0.1  atmosphere  exhibit  excellent  agree¬ 
ment  with  the  JANAF  phase  equilibrium  predictions.  At  the  1.0  atmosphere  condition,  only  the  highest 
temperature  data  point  is  included  in  Figure  10  since  the  other  data  are  obviously  erroneous  due  to 
radiation  blockage.  The  highest  temperature  data  point  is  about  230°R  below  the  JANAF  vapor  phase 
equilibrium  temperature  at  the  1.0  atmosphere  condition,  whereas  it  is  above  both  the  Dolton,  et  al. 
and  Kratsch,  et  al.  values  by  about  430°R  and  250°R  respectively.  Interestingly  it  is  only  about 
70°R  below  the  JANAF/Strauss-Thiele  phase  equilibrium  prediction.  It  is  likely  that  this  surface 
temperature  measurement  is  also  low  since  although  the  plume  was  visually  reduced  by  the  cold  argon 
flow  and  the  N02  concentration  was  obviously  reduced  both  were  still  discernible.  In  addition  the 
supersonic  flow  surface  temperature  data  are  about  30°R  above  the  peak  subsonic  flow  surface  temper 
ature  while  the  B1  values  are  about  a  factor  of  4  lower.  Although  it  is  difficult  to  quantify,  all 
available  data  suggest  cnat  all  of  the  surface  temperature  data  from  the  subsonic/atmospheric  tests 
are  low.  Considering  this,  the  peak  surface  temperature  value  at  one  atmosphere  tends  to  verify  the 
JANAF  thermochemistry  model  since  the  measured  value  is  probably  lower  then  the  true  surface  temper¬ 
ature.  Certainly,  when  taken  as  a  complete  set  the  equilibrium  sublimation  data  support  the  nominal 
JANAF  thermochemical  data. 

4.2  DATA  COMPARISON  WITH  EQUILIBRIUM  AND  RATE  CONTROLLEO  M00ELS 

The  analyses  in  this  section  compare  the  supersonic  flow  (high  mass  transfer  rate)  ATJ-S 
graphite  ablation  data  presented  in  Table  10  with  equilibrium  and  kinetically  controlled  sublimation 
predictions.  The  ablation  data  generated  under  the  high  convective  conditions  are  compared  with  the 
appropriate  JANAF  equilibrium  diffusion  controlled  B1  curve  in  Figure  20.  The  theoretical  B1  curves 
included  in  Figure  20  span  the  pressure  range  0.2  to  1.0  atmospheres  as  do  the  data.  The  impact 
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Figure  20.  Comparison  of  ATJ-S  ablation  data  from  the 
predictions. 


pressures  at  the  supersonic  flow  conditions  for  the  low  (~  0.2  atmosphere)  and  moderate  (~  1  atmo¬ 
sphere)  pressure  tests  are  in  the  ranges  0.2  to  0.28  atmospheres  and  0.8  to  0.9  stmospneres,  respec¬ 
tively.  Therefore  theoretical  8/  curves  for  these  pressures  are  included  in  Figure  20  to  afford  a 
more  accurate  comparison. 

Comparison  of  the  high  convective  B'  data  with  the  JANAF  equilibrium  theory  shows  that  the 
experimental  data  are  consistently  at  slightly  higher  temperatures  than  the  theoretical  B'  curves. 
Since  the  agreement  between  the  equilibrium  B'  data  (i.e.,  B'  »  1)  and  the  JANAF  phase  equilibrium 
is  good  the  ablation  data  comparison  in  Figure  20  shows  a  slight,  but  consistent  shift  to  higher 
ablation  temperatures  is  occurring  with  increased  heat  and  mass  transfer  rates. 


Figure  21  shows  plots  of  the  temperature  shift  (defined  as  Ts  -  TJANAF)  as  a  function  of 
i'(,u0  r'm  for  tno  ,ow  Prer>s,lro  tests.  The  results  in  Figure  21  show  that,  while  the  temperature 
data  al.  a  given  mass  transfer  coefficient  condition  display  only  about  a  2  percent  scatter,  the  \T 
associated  with  this  scatter  is  similar  to  the  AT  anticipated  due  to  sublimation  kinetics.  This  re¬ 
sult  makes  firm  conclusions  regarding  kinetic  effects  difficult.  However,  considering  the  low  pres¬ 
sure  data  in  Figure  21  in  a  statistical  fashion,  it  is  noted  that  there  is  an  average  temperature 
shift  of  about  +65°R  with  increasing  mass  transfer  coefficient.  This  average  shift  is  about  30  per¬ 
cent  of  the  predicted  shift  based  on  the  vaporization  coefficients  recommended  by  Dolton,  et  al. 
(Reference  22),  based  on  the  work  of  Zavitsanos  (Reference  7)  fc  -  the  C1  -  C5  vapor  species. 

Figure  22  shows  the  surface  temperature  shifts  at  the  high  convective  flow  conditions  mea¬ 
sured  from  the  average  of  the  experimental  "equilibrium"  surface  temperatures  at  the  low  pressure 
test  conditions.  Included  in  this  figure  are  the  predicted  temperature  shifts  assuming  unity  vapor¬ 
ization  coefficients  for  each  of  the  carbon  vapor  species  Cj  -  Cg.  It  is  interesting  to  note  that 
the  bulk  of  the  data  lie  between  the  two  sets  of  predictions.  Lundell  and  Dickey  (Reference  47) 
recently  completed  a  graphite  sublimation  test  program  using  high  intensity  (25  K  watt)  C02  laser 
radiation.  To  interpret  their  ablation  results  they  assumed  unity  vaporization  coefficients  for 
the  vapor  species  Cj  -  Cfj  which  resulted  in  good  agreement  with  predictions  using  the  JANAF  data. 

The  results  in  Figure  22  show  the  unity  vaporization  coefficient  kinetic  sublimation  predictions  to 
be  somewhat  below  the  experimental  data,  but  only  about  10°R  -  15"R  below  the  average  shift  of  about 
65"R  at  a  mass  transfer  rate  of  0.1  lbm/ft7  sec. 


The  data  plotted  in  Figure  23  from  the  moderate  pressure  tests  also  exhibit  a  positive  tem¬ 
perature  shift  with  increasing  values  of  pgue  C  .  Interestingly  the  magnitude  of  the  shift  at 
these  conditions  (P  =  0.8  to  1.0  atmospheres)  is  in  the  range  50°R  to  150°R,  somewhat  lower  than 
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Figure  21 . 


Measured  temperature  shift  with  peueCm  for  ATJ-S  ablation  tests 
in  the  pressure  range  of  0.1  to  0.3  atmospheres.  (Referenced 
to  theoretical  equilibrium  temperatures). 
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Figure  23.  Measured  temperature  shift  with  peup('m  for  ATJ-S  ablation 
tests  in  the  pressure  range  0.8  to  i.U  atmospheres. 
(Referenced  to  theoretical  equilibrium  temperatures). 
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the  shifts  at  the  low  pressure  test  condition.  This  trend  of  a  decreasing  temperature  shift  with 
P(up  C(|  with  increasing  pressure  is  in  agreement  with  the  Knudsen-Langmuir  predictions 

Recause  the  experiiiental  “equilibrium"  temperatures  are  felt  to  be  low  due  to  radiation 
blockage  problems  at  the  moderate  pressure  test  conditions  a  temperature  shift  referenced  to  an 
average  experimental  "equilibrium"  temperature  is  not  possible.  Thus,  with  the  limited  data  at 
this  conditioa  one  can  only  observe  that  the  experimental  shift  (referenced  to  a  theoretical 
equilibrium  value)  exhibits  good  agreement  v  th  the  predicted  temperature  shift  due  to  sublimatiom 
kinetics  (assuming  nominal  vaporization  coefficients  for  the  species  -  C^).  It  is  important  to 
note,  however,  that  the  temperature  shifts  at  this  moderate  pressure  condition  are  within  the  +2 
percent  uncertainty  of  the  temperature  data.  Thus,  it  is  not  practical  to  draw  any  firm  conclusions 


SECTION  5 

ABLATION  DATA  COMPARISON  WITH  PREVIOUSLY  PUBLISHED  RESULTS 

This  section  compares  the  ATJ-S  ablation  data  from  this  test  series  with  the  ablation  data 
from  previous  investigations.  The  results  of  two  ablation  test  programs,  both  conducted  in  the 
Ames  AEHS  facility  are  reviewed  and  ablation  data  comparisons  made.  Data  from  Wakefield  and  Petersons 
(Reference  48)  tests  of  ATJ  graphite  in  the  AEHS  facility  are  reviewed  and  compared  with  the  low 
pressure  supersonic  flow  ablation  data  from  this  series  in  Section  5.1.  Oavy  and  Bar  Nun  (Ref¬ 
erence  38)  also  studied  graphite  ablation  under  radiation  only  test  conditions  in  the  AEHS  facility. 
Their  test  results  are  reviewed  and  ablation  data  comparisons  made  in  Section  5.2.  In  addition, 

AT.l  ablation  data  from  the  nominally  one  atmosphere  ablation  tests  of  Lundell  and  Oickey  (Reference 
24)  which  were  conducted  in  the  NASA-Ames  Heat  Transfer  Tunnel  are  reviewed  and  compared  with  data 
selected  from  this  program  in  Section  5.3. 

5.1  WAKEFIELD  AND  PETERSON'S  AEHS  ABLATION  TEST  SERIES 

Wakefield  and  Peterson's  ablation  tests  of  hemi spherically  tipped  hollow  ATJ  graphite  ablation 
models  were  similar  to  this  test  series  with  the  following  exceptions: 

1.  Wakefield  and  Peterson's  ablation  tests  were  limited  to  the  0.06  to  0.33  atmosphere 
supersonic  flow  regime. 

2.  The  Linde  arc  heater  was  used  in  the  AEHS  facility  instead  of  the  constrictor  heator  used 
on  this  program. 

3.  The  optical  pyrometer  used  by  Wakefield  and  Peterson  was  a  monochromatic  device  (Thermodot 
T0-9)  with  a  narrow  bandpass  filter  centered  at  0.8p. 

The  difference  between  the  optical  pyrometers  used  on  these  two  test  series  provides  an  inter¬ 
esting  check  on  the  accuracy  of  the  pyrometer  data  from  both  test  series.  The  response  time  of  the 
TD-9  pyrometer  is  reported  to  be  about  100  msec  for  a  30  percent  full  scale  deflection.  Thus  one 
would  not  expect  the  pyrometer  signal  to  respond  to  the  3  msec  interval  during  which  the  incident 
radiation  is  blocked.  (Interestingly,  Wakefield  and  Peterson  report  that  the  component  of  reflected 
radiation  is  small  as  one  would  expect  with  the  TD-9  response.)  Results  reported  in  Section  3.2.3 


76 


showed  the  Thermogaqe  pyrometer  recorded  about  a  200°R  temperature  drop  during  the  3  msec  the  incident 
radiation  was  blocked.  One  must  consider,  however  that  the  optical  sensitivity  of  the  Thermoaage 
covers  a  broader  spectral  range  than  that  of  the  TD-9,  thus  its  sensitivity  to  the  reflected  radiation 
is  greater  (note  the  spectral  distribution  of  the  AEHS  radiation  source  in  Figure  8).  In  addition 
to  these  differences  in  the  optical  pyrometers,  Wakefield  and  Peterson  corrected  the  measured  bright¬ 
ness  temperatures  by  assuming  a  surface  emissivity  of  0.9  This  corresponds  to  an  "actual"  surface 
temperature  about  130"R  above  the  brightness  (.  =  1.0)  temperature  for  these  conditions.  Wakefield 
and  Peterson  report  a  3  percent  uncertainty  in  the  surface  temperature  data,  which  corresponds  to 
about  a  200  R  uncertainty.  A  minimal  amount  of  mass  loss  in  the  condensed  phase  was  observed  by 
Wakefield  and  Peterson,  based  on  their  interrogation  of  film  data  and  observation  of  the  post-test 
models. 

The  comparison  between  Wakefield  and  Peterson's  results  and  the  0.22  to  0.28  atmosphere  super¬ 
sonic  flow  results  from  this  test  series  compared  in  Fiqure  24  shows  excellent  agreement.  The  surface 
temperature  data  in  Fiqure  24  exhibits  agreement  within  <200°R  which  is  within  the  basic  data  uncertainty. 
Wakefield  and  Petersen  compared  their  test  data  with  ablation  predictions  based  on  JANAF  thermochemistry 
for  the  carbon  vapor  species  and  found  excellent  agreement.  This  correlation,  however  did  not  conclu¬ 
sively  establish  the  JANAF  t.hertnochemi stry  for  the  carbon  vapor  species  as  being  correct,  since  the 
tests  were  conducted  under  relatively  high  heat  and  mass  transfer  rate  conditions  which  could  imply 
significant  sublimation  kinetic  effects.  The  excellent  corroboration  between  their  ATJ  ablation  data 
and  the  ATJ-S  data  from  this  test  series  further  substantiates  the  accuracy  of  the  reduced  data 
previously  reviewed  in  Section  3.4.3. 

5.2  0AVY  AN0  BAR  NUN'S  AEHS  ABLATION  TEST  SERIES 

This  ablation  test  series  in  the  AEHS  facility  tested  ATJ  graphite  models,  in  the  form  of 
solid  0.4  inch  diameter  right  circular  cylinders  1.5  inch  long.  The  ablation  tests  were  radiation 
only  tests  in  a  cold  argon  environment.  The  test  chamber  environment  was  prepared  for  each  test  by 
pumping  down  to  pressures  of  1-10,.  for  10-20  minutes,  purging  the  chamber  with  argon,  and  then  filling 
with  room  temperature  argon  to  the  desired  test  pressure. 

Davy  and  Bar  Nun  used  the  same  Thermodot  T0-9  FH  optical  pyrometer  as  Wakefield  and  Peterson. 

Ouring  the  test  series  they  made  both  front  and  sideview  readings  of  the  ablating  models  at  various 
times  at  comparable  test  conditions  and  observed  no  essential  differences  in  the  ablation  temper¬ 
ature.  They  also  reported  that  after  the  initial  heating  transient,  the  pyrometer  signal  reached  a 
steady  state  condition  in  which  the  chopper  modulation  was  only  barely  perceptible.  This  result  is 
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Wakefield  &  Peterson's  data  Data  from  this  series 


in  qood  agreement  with  that  reported  by  Wakefield  and  Peterson  as  would  be  expected.  At  pressures 
above  0.1  atmosphere  slight  temperature  corrections  were  made  due  to  soot  blockage  of  the  incident 
and  emitted  surface  radiation.  In  the  pressure  range  0.1  to  0.3  atmospheres  these  soot  corrections 
to  the  measured  temperatures  were  less  than  90°R.  At  the  1.0  atmosphere  test  condition  the  soot 
correction  to  the  measured  surface  temperature  was  about  500°R. 

The  radiation  only  data  from  Davy  and  Bar  Nun's  ATJ  ablation  tests  are  compared  with  the  low  sub¬ 
sonic  flow  test  data  from  this  series  in  Figure  35.  Davy  and  Bar  Nun's  data  in  Figure  25  are  plotted 
using  the  soot  corrected  surface  temperatures  and  the  assumption  that  the  equilibrium  carbon  vapor 
pressure  is  equivalent  to  the  static  chamber  pressure.  Despite  Uavy  and  Bar  Nun's  attempt  to  correct 
the  measured  temperatures  (due  to  soot  absorption)  their  data  are  consistently  about  400°R  below 
the  data  from  this  series.  The  high  enthalpy  low  subsonic  flow  employed  for  these  tests  removed 
any  potential  soot  film  at  the  0.2  to  0.3  atmosphere  test  conditions.  At  the  1.0  atmosphere  test 
condition  (cold  subsonic  argon  flow)  some  blockage  was  felt  to  exist  durinq  these  tests  but  no  attempt 
was  made  to  quantify  this  effect.  Despite  the  fact  that  the  measured  surface  temperature  from  this 
test  series  at  the  1.0  atmosphere  test  condition  is  considered  low,  it  is  still  about  400°R  above 
the  highest  value  reported  by  Davy  and  Bar  Nun.  One  can  only  conclude  from  this  comparison  that 
the  soot  absorption  which  occurred  due  to  the  cold  argon  environment  was  underestimated.  In  addition, 
the  blunt,  solid  ATJ  graphite  ablation  models  tested  by  Davy  and  Bar  Nun  provided  a  significantly 
higher  capacitance  than  the  models  tested  on  this  series.  This  difference  would  cause  lower  steady 
state  ablation  temperatures,  which  is  substantiated  by  the  fact  that  the  mass  loss  rates  reported  by 
Davy  and  Bar  Nun  are  about  an  order  of  magnitude  lower  than  the  data  from  this  series. 

Based  on  this  ablation  data  comparison  and  differences  in  the  two  test  techniques  one  can 
identify  several  reasons  why  the  ablation  temperatures  measured  by  Davy  and  Bar  Nun  would  be  low 
Therefore,  the  conclusions  regarding  the  carbon  vapor  thermochemistry  drawn  by  these  researchers  are 
highly  suspect.  There  still  remains  the  question  regarding  the  degree  of  incident  and  emitted 
radiation  blockage  which  occurred  in  this  test  series  at  the  one  atmosphere  test  condition.  Sugges¬ 
tions  as  to  possible  test  techniques  available  to  further  study  this  problem  are  reviewed  in  Section  8. 

5.3  LUNDELl  AND  DICKEY'S  ATJ  ABLATION  TEST  SERIES 

The  ATJ  ablation  test  series  conducted  by  Lundell  and  Dickey  (Reference  24)  has  become  a 
classical  reference  within  the  carbon  ablation  community.  This  relatively  extensive  ablation  test 
series  covered  the  pressure  range  of  0.3  to  4.4  atmospheres  and  was  conducted  in  the  NASA  Ames  Heat 
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Transfer  Tunnel.  The  intent  of  their  test  program  was  the  generation  of  high  quality  ATJ  ablation 
data  in  the  sublimation  regime.  A  further  objective  of  their  program  was  to  determine  if  significant 
ablation  mechanisms  other  than  the  classical  thermochemical  phenomena  exist  ( e . g . ,  particulate  mass 
loss). 


The  Aines  Heat  Transfer  Tunnel  is  a  conventional  arc-driven  wind  tunnel.  For  these  tests  a  Linde 
arc  heater  (Model  N4001  )  was  used.  The  flow  over  the  models  varied  from  Mach  3.1  to  3.8.  During 
a  run  as  many  as  18  separate  support  arms  are  available  for  insertion  into  the  flow.  Instrumentation 
inserted  into  the  test  stream  during  each  run  included  calorimeters  (both  slug  type  and  steady  state) 
and  pressure  probes.  Instrumentation  external  to  the  test  chamber  included  optical  pyrometers  and  a 
motion  picture  camera.  Two  pyrometers  were  used  so  that  the  monochromatic  brightness  temperature  would 
be  measured  at  two  different  wavelengths.  The  pyrometers  used  were  a  Thermodot  Model  TD-9  with  a 
bandpass  at  0.8p  and  a  Pyro  650  with  a  bandpass  filter  at  0.65p.  Both  pyrometers  were  calibrated  to 
6 1 00 " R  through  use  of  a  conventional  ciisappearing-filament  optical  pyrometer. 


Lundell  and  Dickey's  ATJ  ablation  data  at  a  nominal  surface  pressure  of  1  atmosphere  are  shown 
in  Figure  26.  The  raw  ablation  data  were  transformed  into  B'  values  through  use  of  the  following 
equations.  The  nondimensional  mass  transfer  parameter  B'  is  defined  as,  B'  =  ^DeUe^H  w^ere: 


-e"eCH„  '  “-113  VV«„ 


Reference  24 


and 


1 . 4Bo 


CH/CH0  =  (e1 • "Bo  -  1) 

The  nose  radius  (Ri  used  in  the  expression  for  the  convective  heat  transfer  coefficient  (p  u  Cu  ) 
n  e  e  no 

was  the  effective  value  which  Lundell  evaluated  and  reported  for  each  ablation  model  (based  on  the 
initial  and  final  nose  radius).  Included  in  Figure  26  is  the  predicted  JANAF  equilibrium  and  kinetic 
sublimation  B'  curves.  The  agreement  between  the  equilibrium  prediction  and  the  data  is  excellent. 
Also  included  in  Figure  26  are  the  1  atmosphere  ablation  data  from  this  test  series.  These  data  are 
at  B  values  nearly  a  factor  of  two  above  the  highest  Lundell  and  Dickey  value.  This  data  comparison 
shows  the  measured  ablation  temperatures  on  this  test  series  to  be  about  300°R  above  the  values  of 
Lundell  and  Dickey.  The  ablation  data  from  this  series  exhibit  good  agreement  with  the  kinetic  sub¬ 
limation  prediction  in  Figure  26. 

The  convective  flow  environments  over  the  models  in  both  tests  are  extremely  similar  as  shown 
by  the  data  in  Table  1 1 . 
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Temperature  *  " 

of  ablation  data  from  this  series  with  the  data 
and  Dickey. 


TABLE  11.  COMPARISON  OF  1  ATMOSPHERE 
SUPERSONIC  F10W  CONVECTIVE 
ENVIRONMENTAL  PARAMETERS. 


This  Series 

Lundell  S  Dickey's 

GdS 

air 

ai  r 

18,000 

-17.000 

Mach  no. 

-2.6 

-3.5 

Model  size  (in.) 

D  =  0.5 

D  =  0.6  to  1.2 

The  principal  difference  between  the  two  test  series  results  from  the  high  radiant  flux  in  .he  AtllS 
facility.  The  temperature  data  in  Figure  26  from  this  test  series  have  been  corrected  for  the  reflected 
component  of  the  incident  radiation.  Similar  ablation  data  from  the  low  pressure  (-0.1  atm.)  super¬ 
sonic  flow  condition  were  compared  with  the  results  of  Wakefield  and  Peterson  in  Section  5.1  and 
shown  to  exhibit  excellent  agreement.  Thus  the  data  quality  from  this  series  is  substantiated  by 
these  direct  comparisons  with  other  ablation  data. 

Lundell  and  Dickey  did  report  a  significant  amount  of  particulate  mass  loss  from  their  blunt 
cylindrical  models.  Althouqh  the  cause  of  this  particulate  mass  flux  is  not  understood,  the  effect 
of  particulate  mass  loss  on  the  surface  energy  balance  is  to  lower  the  ablation  temperature.  There¬ 
fore,  the  most  probable  explanation  for  the  300°R  surface  temperature  difference  between  these  two  sets 
of  ablation  data  is  the  particulate  mass  loss  which  Lundell  and  Dickey  nhserved  from  their  models 
in  the  Ames  Heat  Transfer  Tunnel.  The  data  from  this  test  series  suggest  negligible  particulate 
mass  loss  occurred  which  corroborates  well  with  the  results  of  Wakefield  and  Peterson.  The  surfaces 
of  the  ablation  models  were  all  smooth  indicative  of  no  gross  spallation  and  review  of  the  hign  speed 
motion  picture  data  showed  virtually  no  particulate  mass  loss.  The  only  particulate  mass  flux  evident 
from  the  high  speed  film  data  occurred  during  Run  81  which  was  unique  since  cold  argon  was  flowing 
over  the  radiatively  heated  ablation  model.  Close  examination  of  this  post-test  model  showed  a  rougn 
irregular  surface  around  the  periphery  of  the  front  ablative  surface  where  the  particulate  mass  loss 
obviously  occurred.  In  this  rather  singular  case  the  severe  thermal  gradients  in  the  region  of  the 
interface  between  the  front  ablating  surface  and  the  cooler  cylindrical  afterbody  resulted  in  material 
spallation.  A  similar  nnenomenon,  but  less  severe  that  this  in  Lundell  and  Dickey's  ablation  tests 
may  have  caused  the  rather  significant  particulate  mass  flux  which  they  observed.  Thus,  the  lack  of 
agreement  between  these  two  sets  of  data  is  best  rationalized  by  the  particulate  mass  flux  which  was 


SECTION  6 


IMPLICATIONS  OF  SUBLIMATION  KINETICS  ON  HIGH  PRESSURE  ABLATION  RESPONSE 


The  AT J -S  ablation  data  analyses  in  Section  4  show  that  the  kinetic  sublimation  predictions 
(using  the  JANAF  thermochemistry  data  and  the  vaporization  coefficients  of  Zavitsanos  (Reference  7) 
extended  by  Dolton,  et  al.  (Reference  22))  slightly  overpredict  the  measured  temperature  shifts  from 
the  theoretical  equilibrium  ablation  temperatures.  These  analyses  were  restricted  to  the  low  and 
moderate  pressure  regime  since  the  ablation  data  were  restricted  to  this  regime. 


The  intent  of  this  section  is  to  extrapolate  these  low  and  moderate  pressure  ablation  results 
to  the  high  pressure  regime  of  interest.  The  kinetically  controlled  sublimation  predictions  are 
made  using  the  JANAF  thennochemi cal  properties,  the  Knudsen-Langmuir  kinetic  model,  and  the  vaporiza¬ 
tion  coefficients  of  Zavitsanos  (Reference  7)  extended  to  C4  and  C5  by  Dolton,  et  al.,  (Reference  22). 
The  thennochemi  cal  1)'  curves  presented  in  this  section  were  predicted  with  the  modified  version  of 
the  Aerotherm  Chemical  Eauilibrium  (ACE)  computer  code  (Reference  35).  The  assessment  of  the  high 
pressure  sublimation  kinetic  effects  are  made  for  typical  AFFDL  50  MW  ablation  environmental  condi¬ 
tions  and  a  typical  flight  environment,  in  addition,  analyses  are  summarized  which  illustrate  the 
sensitivity  of  kinetic  sublimation  predictions  to  the  state  of  the  carbon  vapor  species  (i.e.,  whether 
frozen  or  equilibrium).  In  total,  these  analyses  illustrate  the  potential  effect  of  carbon  sublima¬ 
tion  uncertainties  in  the  high  pressure  regime  (i.e.,  p  >  10  atm)  as  extrapolated  from  the  low  pres¬ 
sure  ablation  data  from  this  program's  AEHS  facility  test  series. 


The  generalized  kinetic  sublimation  predictions  are  reviewed  in  Section  6.1  where  comparisons 
are  made  with  the  equilibrium  predictions  and  the  sensitivity  of  kinetic  predictions  to  the  carbon 
vapor  state  are  studied.  Results  of  analyses  showing  the  effect  of  uncertainties  in  the  carbon 
vapor  state  are  reviewed  in  substantially  more  detail  than  the  cursory  review  presented  in  Section 
2.3.3.  Section  6.2  then  compares  the  effect  of  sublimation  kinetics  on  the  predicted  ablation  re¬ 
sponse  of  selected  high  pressure  50  MW  ablation  data  and  a  typical  reentry  condition. 
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6.1  GENERALIZED  KINETIC  SUBLIMATION  ABLATION  PREDICTIONS 

Hip  kinetic  sublimation  predictions  presented  in  this  section  are  restricted  to  moderate  and 
hi.jh  pressures  (i.e.,  1  -  300  atm)  since  over  90  percent  of  nosetip  recession  generally  occurs  in 
this  pressure  regime.  These  predictions  are  designed  to  answer  the  question:  What  is  the  most 
probable  effect  of  carbon  sublimation  kinetics  on  the  predic‘ed  ablation  response  for  current  envi¬ 
ronmental  conditions  of  interest?" 

Results  of  these  calculations  in  the  form  of  B 7 temperature  plots  are  shown  in  Figure  27, 
Pressures  from  1  to  300  atmospheres  and  mass  transfer  coefficients  covering  the  range  of  0.1  to  10 
lbm/ft2sec  were  assumed  to  generate  these  results.  The  carbon  vapor  species  were  assumed  to  be  either 
frozen  (as  was  assumed  in  Section  4)  or  to  equilibriate  in  the  vapor  state.  The  frozen  vapor  predic 
tions  are  used  for  comparison,  however  since  this  assumption  maximizes  the  shift  from  the  equilibrium 
diffusion  controlled  ablation  state. 

It  is  clear  from  these  results  that  at  high  pressure  ablation  conditions  of  interest  (i.e.. 

Pl?  •  50  atmospheres)  the  uncertainty  in  the  surface  temperature  resulting  rrom  an  uncertainty  in  the 
sublimation  response  of  carbon  is  probablv  below  100"R  (based  on  an  extrapolation  of  the  data  analyses 
in  Section  4  and  the  high  pressure  predictions  shown  in  Figure  27). 

Not  only  does  the  kinetic  sublimation  response  of  graphitic  materials  present  an  uncertainty 
in  ablation  predictions  (although  the  results  in  Section  4  put  a  bound  on  the  kinetic  effects)  but 
the  state  of  the  carbon  vapor  species  at  the  surface  continues  to  be  uncertain  as  discussed  in  Sec¬ 
tion  2.  Early  kinetic  sublimation  models  assumed  the  carbon  vapor  species  to  be  frozen.  More  re¬ 
cent  kinetic  sublimation  models  have  assumed  the  carbon  vapor  to  be  in  equilibrium.  A  direct  compar¬ 
ison  of  these  two  assumptions  has  been  made  to  assess  the  effect  of  the  carbon  vapor  chemical  state 
on  the  predicted  ablation  response.  These  results  are  shown  in  Figure  27.  The  comparison  is  made 
for  a  mass  transfer  coefficient  of  10  lbm/ft2sec  since  sublimation  kinetic  effects  are  greatest  at 
high  mass  transfer  conditions. 

The  significance  of  the  assumed  carbon  vapor  state  increases  with  temperature  and  pressure. 
Little  difference  exists  between  the  equilibrium  and  frozen  carbon  vapor  solutions  at  1  atmosphere, 
whereas  at  300  atmospheres  the  equilibrium  vapor  case  corresponds  more  closely  to  the  equilibrium/ 
diffusion  controlled  solution  than  the  kinetic  sublimation/frozen  vapor  solution.  This  results  for 
the  following  reasons.  When  the  vapor  state  is  frozen  each  kinetic  reaction  controls  a  portion  of 
the  total  sublimation  response.  However,  when  the  vapor  species  can  equilibrate  the  sublimation  re¬ 
sponse  becomes  controlled  by  only  the  fastest  kinetic  reaction.  This  shift  is  maximized  at  high  pres¬ 
sures  because  of  the  increased  significance  of  the  higher  order  vapor  species  (e.g.,  Cg)  which  has  an 
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27.  Concluded. 


extremely  small  vaporization  coefficient.  Thus,  when  the  vapor  is  frozen,  kinetic  sublimation  of  tg 
is  (, really  restricted.  Conversely,  when  the  vapor  can  equilibrate  the  partial  pressure  potential 
ol  the  faster  sublimation  reactions  is  increased  due  to  the  gas  phase  reactions  which  increase  the 
concentration  of  those  vapor  species  restricted  by  their  slow  sublimation  rates  (i.e.,  C3  and  C&) 
while  reducing  the  concentration  of  the  faster  subliming  species  (i.e.,  C3). 

This  is  readily  shown  by  the  predicted  ablation  data  summarized  in  Table  12.  At  the  300  atmo¬ 
sphere  condition  the  kinetic  sublimation  rates  of  C?  and  C3  change  by  factors  of  5  and  1/3  respec¬ 
tively  when  shifting  from  an  assumed  frozen  to  equilibrium  vapor  at  the  surface.  Thus,  the  equilibrium 
vapor  assumption  causes  the  C?  concentration  to  be  reduced  while  the  C3  and  C5  concentrations  are  sub¬ 
stantially  increased  (particularly  at  higher  temperatures).  For  this  reason,  kinetic  sublimation  pre¬ 
dictions  based  on  an  assumed  equilibrium  carbon  vapor  at  the  surface  are  shifted  back  toward  the  equi¬ 
librium/diffusion  controlled  predictions  as  the  temperature  increases  (i.e.,  higher  pressures). 

Relating  this  result  to  the  ablation  data  generated  in  this  program  in  a  precise  quantitative 
manner  is  not  possible.  However,  the  data  analyses  in  Section  4  show  a  consistent  shift  to  tempera¬ 
tures  above  the  equilibrium  diffusion  controlled  temperature  predictions  under  high  heat  and  mass 
transfer  conditions.  The  data  from  this  program  however  are  insufficient  to  back  out  either  vapor¬ 
ization  coefficients  or  to  assess  whether  the  carbon  vapor  species  are  frozen  or  chemically  equili¬ 
brate.  The  data  suggests  an  upper  bound  of  these  chemical  kinetic  effects  and  this  is  evaluated  in 
the  next  section. 

6.2  KINETIC  SUBLIMATION  PREDICTIONS  APPLIED  TO  50  MW  ABLATION  DATA  AND  FLIGHT 

To  quantitatively  assess  the  effects  of  sublimation  kinetics  on  the  predicted  ablation  response 
ol  graphitic  materials,  kinetic  predictions  were  compared  with  cquilihrium/diffusion  controlled  pre¬ 
dictions  of  recent  50  MW  ATJ-S  arid  C/C  tests.  The  kinetic,  predictions  assumed  a  frozen  carbon  vapor 
in  order  to  maximize  kinetic  effects.  Only  the  high  pressure  turhulent  tests  were  analyzed  since  the 
ahlation  response  and  temperature  measurements  are  both  definitely  steady  state  and  easily  related 
to  the  ablation  model  due  to  the  turbulent  biconic  shape  which  develops.  Rough  wall  heating,  based 
on  a  micro  roughness  height  of  1  mil  were  assumed  for  these  analyses.  The  correction  factor  (Kr)  to 
the  smooth  wall  heat  and  mass  transfer  coefficient  was  based  on  the  correlation  of  Powars  (Reference  49). 

The  kinetic  and  equilibrium/diffusion  controlled  predictions  are  compared  in  Figure  2B.  Both 
suh I imation  modeling  techniques  are  seen  to  predict  the  measured  ablation  temperature  within  a  200  R 
(i.e.,  '2.5  percent)  uncertainty  band.  Kinetically  controlled  sublimation  predicts  ablation  tempera 
tures  about  45”R  (0.5  percent)  above  the  diffusion  controlled  values.  The  decrease  in  the  predicted 
recession  rate  with  sublimation  kinetics  is  less  than  2  percent  for  the  test  conditions. 
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TABLE  12.  SENSITIVITY  OF  CARBON  VAPOR  SPECIES  CONCENTRATIONS  AND  SUBLIMA¬ 
TION  RATES  WITH  THE  ASSUMED  CARBON  VAPOR  STATE  ANn  WITH  PRESSURE 


Pressure  -  1  atm 
T  =  7605°R 


State  of 

Carbon  Vapor 

Frozen 

Equil ibr ium 

Bc 

C1 

.02722 

.0186 

Bc 

c2 

.11756 

.13784 

Bc 

c3 

.04481 

.04408 

Bc 

C4 

.00801 

.00877 

Bc 

c5 

.00017 

.00017 

ks 

.1978 

.2136 

p  • 

“total 

1.0 

1.018 

Species 

Mole  Fraction 

State  of  Carbon  Vapor 

Frozen 

Equilibrium 

C1 

.3387-1 

.6204-1 

C2 

.7205-1 

.2929-1 

C3 

.1856-1 

.4663-1 

C4 

.2515-2 

.2172-3 

C5 

.6771-4 

.1502-3 

CN 

.2259 

.2205 

c2n 

.2694 

.2669 

C2N2 

.4516-2 

.4432-2 

C4N2 

.1280-2 

.1256-2 

CO 

.2181 

.2153 

Pressure  -  300  atm 
T  =  8993#R 


State  of 

Carbon  Vapor 

Mole 

Fraction 

Frozen 

Equilibrium 

State  of 

Carbon  Vapor 

B ' 

.0112 

.0186 

Frozen 

Equilibrium 

C1 

B' 

.0553 

.2689 

C1 

.1507-1 

.1490-1 

C2 

C2 

.3707-1 

.3544-1 

Bc 

.2510 

.0896 

C3 

.1120 

.1331 

B' 

.0145 

.0508 

C4 

.4871-2 

.4489-2 

c4 

C5 

.2139-2 

.1032-1 

B' 

.0079 

.0013 

c5 

CN 

.1035 

.1012 

5 

c2n 

.2046 

.2002 

Vb' 

.3399 

.4292 

3“T  cj 

c2n2 

.2655-1 

.2542-1 

* 

C4N2 

.1990-1 

.1905-1 

p  • 

“total 

1.0 

1.094 

CO 

.2342 

.2259 

C/C  data/predictions 
ATJ-S  data/predictions 


Figure  28.  Effect  of  sublimation  kinetics  on 
predicted  graphite  ablation  temper¬ 
atures  in  the  50  MW  arc. 


SECTION  7 


CONCLUSIONS 

Conclusions  regarding  the  thermochemical  sublimation  response  of  ATJ-S  graphite  derived 

from  high  temperature  ablation  tests  under  low  and  high  mass  transfer  environments  in  the  NASA 

Ames  AEHS  facility  are  summarized  below. 

1.  The  variable  velocity  flow  (subsonic  to  supersonic)  test  procedure  designed  for  the  NASA  Ames 
AEHS  facility  proved  to  be  a  viable  test  technique  for  studying  the  detailed  thermochemistry 
of  graphite  sublimation  in  the  low  to  moderate  pressure  regime. 

?.  The  high  IV,  low  I'pU  C^  equilibrium  data  generated  on  this  program  exhibit  very  qood  aqreement 

with  the  carbon  vapor  pressure  curve  predicted  by  the  JANAF  thermochemical  data  for  carbon 

vapor  species  C,  through  C^. 

3.  The  ablation  data  generated  on  this  program  displayed,  on  the  average,  an  increasing  shift  in 

surface  temperature  above  the  theoretical  (JANAF)  equilibrium  value  with  increasing  convective 
mass  transfer  coefficient.  This  shift  is  in  the  direction  to  substantiate  that  the  graphite 

ablation  rate  is  being  partially  controlled  by  sublimation  kinetics. 

4.  The  measured  temperature  shifts  with  increasing  mass  transfer  coefficient  were  relatively 

small.  At  the  low  pressure  condition  (0.22  to  0.28  atm.),  the  averaoe  shift  in  T  .  - 

measured 

^ JANAF  was  about  150°R  as  PeueCM  increased  from  near  zero  to  0.08  -  0.11  lbm/ft2sec.  The 
shift  predicted  using  the  kinetic  sublimation  coefficients  reported  by  Zavitsanos  as  extended  by 
Dolton,  et  al.,  for  these  conditions  is  ~200°R.  Thus,  it  may  be  concluded  that  graphite  sub¬ 
limation  kinetics  effects  in  the  moderate  pressure  regime  are  probably  no  greater  than  predicted 
by  Dolton,  et  al . 

5.  While  the  surface  temperature  measurements  were  judged  to  be  quite  accurate  (±2  percent)  the 

temperature  data  scatter  on  an  absolute  scale  was  about  100°R  which  is  similar  to  the  measured 

temperature  shifts.  Thus,  the  data  are  felt  to  be  too  imprecise  to  "back  out'  quantitative 

sublimation  coefficients  or  to  determine  the  thermodynamic  state  of  the  carbon  vapor  near  the 

surface  (i.e.,  equilibrium  or  frozen). 
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6.  Based  on  the  results  of  this  program,  the  most  appropriate  thermochemical  moaei  Tor  grapnite 
nosetip  flight  performance  predictions  is  equilibrium  diffusion  controlled  ablation  using  the 
JANAF  thermochemical  data. 


■  ■- 


SECTION  8 


RECOMMENDATIONS 

The  recommendations  for  further  studies  of  the  thermochemical  ablation  response  of  graphitic 
materials  address  two  areas  of  the  problem.  The  first  addresses  the  more  detailed  thermochemical 
events  which  control  graphite  sublimation  in  air  {much  in  the  same  manner  as  this  program  did). 

The  second  is  a  more  practical  extension  of  the  work  of  Lundell  and  Dickey  (References  24  and  25) 
to  the  higher  pressure  regime  of  interest  for  reentry.  Both  recommendations  are  listed  below. 

1.  The  work  accomplished  on  this  program  has  established  the  AEHS  ablation  test  facility 
as  a  useful  tool  for  studying  graphite  sublimation  under  low  and  high  mass  transfer  en¬ 
vironments.  The  test  techniques  have  been  established  and  thoroughly  checked  out.  The 
data  base  established  on  this  program,  however,  is  not  sufficient  for  meaningful  statis¬ 
tical  analyses.  Additional  ablation  testing  in  the  AEHS  facility  is  recommended  to  pro¬ 
vide  the  data  necessary  for  more  quantitative  conclusions  to  be  drawn. 

2.  During  the  course  of  this  program  a  study  was  conducted  to  assess  the  capability  of 
existing  high  pressure  ablation  test  facilities  with  respect  to  successfully  conducting 
graphite  sublimation  tests  in  the  20  to  40  atmosphere  pressure  range.  Results  of  this 
study  are  documented  in  References  50  and  51.  These  studies  showed  that  all  existing 
facilities  were  inadequate  in  this  regard,  due  mainly  to  enthalpy  limitations  which 
restrict  the  achievable  B‘  levels.  However,  there  are  currently  two  arc  heaters  under  1,3 
construction  which  may  have  the  requisite  pressure  and  enthalpy  capabilities  to  study 
graphite  sublimation  at  these  higher  pressures;  these  are  the  KBC-100  arc  at  AFFDL  and 
the  HEAT  arc  at  AEDC.  It  is  recommended  that  if  the  required  enthalpy  levels  are  demon¬ 
strated  in  either  of  these  facilities  then  moderate  to  high  pressure  graphite  vaporiza¬ 
tion  tests  should  be  conducted.  These  data  would  serve  to  extend  the  results  of  Lundell 
and  Dickey  to  higher  pressures  and  would  have  very  practical  implications  on  graphite  ab¬ 
lation  thermochemistry  in  the  range  of  interest  for  flight  nosetip  ablation  predictions. 
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FOOTNOTES  TO  TABLES  B-l  AND  B-2 
APPENDIX  B 

1.  Initial  system  sheckout  runs:  starting  procedure  verified,  Pgg  rise  calibrated  pumps  required 
defined,  exposure  time  determined,  instrumentation  verified. 

2.  Radiation  source  went  out  during  model  exposure,  restarted  after  models. 

3.  Runs  19  k  21  no  data  available. 

4.  Two  injections,  H.y-Cal  data  repeated. 

6.  Radiation  chopper  not  on  for  this  test.  T^  w/o  RELF  estimated  from  runs  25  &  28. 

6.  Main  valve  closed  data  not  reduced. 

7.  Run  29:  starting  procedure;  run  30:  no  data  available;  run  31:  N2  load  pressure  fell  off  - 
main  supply  closed. 

3.  Exposed  to  combined  convective  and  radiative  flow. 

9.  Run  40:  no  deflection  no  data 

10.  New  slug  c.ilorimeier  installed. 

11.  Medtherm  12  calorimeter  questionable  data. 

12.  Slug  exposed  to  convective  flow,  therefore  not  blacked  for  radiation  after  measurement. 

13.  All  probes  injected  downstream  of  test  station  =1.0". 

14.  Run  47:  no  deflection  on  data  channels. 

15.  Based  on  average  of  after  from  runs  25,  26,  28. 

16  Hg  value  assuming  54  percent  arc  heater  efficiency  (n ) -  No  enthalpy  peaking,  and  fully  filled 

nozzle. 

17.  Initial  Pt2  trace  off  scale.  ?t2  -  A  +  PgOX  @  injection. 

18.  Hg  value  assuming  65  percent  arc  heater  efficiency  (n) ,  no  enthalpy  peaking,  and  fully  filled 
nozzle. 
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I1), 


/II. 


21 


//. 


/I. 

24. 

25. 
25. 

27. 

28. 


29. 


30. 

31. 

32. 


33. 


34. 

35. 
35. 
37. 


i  ■ 


Huif.  55-57  stiirt.  mg  procedure  verified. 

11. Ill  (  I  C.  <|,|  I  VI) 

Hun  59  laser  riol  on,  |>yro  mil  visible  on  (1C,  prool  test.  model  exposed  data  not  reduced. 

Chopper  i‘  f  I  cot.  not  visible  on  0‘graph,  m.ix  temperature  curve  was  ostimatod  from  reflective 
i  omponeii I,  of  similar  runs. 

Huns  63-5fi  starling  procedure  verified. 

Runs  74-75  starting  procedure  verified. 

Reduced  cabin  pressure  for  start  only,  bleed  up  to  1.0  atm  for  test. 

New  slug  calorimeter  installed  for  remainder  of  test. 

N?  only. 

IID  value  assuming  74  percent  arc  heater  effeciency  (n),  no  enthalpy  peaking,  and  fully  filled 
nozzle. 

11^  value  assuming  40  percent  arc  heater  effeciency  (n).  no  enthalpy  pending,  and  fully  filled 
nozzle  and  average  of  voltages  for  this  condition. 

Based  on  average  Pt2  values  from  runs  44  &  45. 

Estimated  from  run  48 

Estimated  from  g  =  50.5  BTU/ft  -sec  on  run  53. 

s  data  from  As/A6)f,  based  on  film  pre-  and  post-profile  measurements.  No  reference  line 
available  for  film  s  reduction. 

Estimated  from  average  ratio  of  s  nom  film/(As/A6)  for  runs  ol  and  63 
Estimated  from  average  of  runs  58,  60-62. 

Estimated  from  previous  run. 

Based  on  average  s  nom/As/Ati  ratio  frew  runs  77  and  79;  76  and  80  not  included  because  no 
s  nom  was  available  and  As/A0  was  significantly  different  than  As/A9)  film. 
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